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sua precisão ao longo do tempo. As células 

de pesagem, devido às suas reduzidas dimensões, 
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Não requerendo as alavancas normalmente usadas 
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As ligações entre as células de pesagem e o aparelho 
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À STANDARD ELECTRICA PROJECTA-SE NO MUNDO 


Eléctrica fabrica numerosos pro- 
dutos, desde os minúsculos tran- 
sistores e circuitos integrados às 


Na vanguarda da industria das Te- 
lecomunicações e Electrônica, a 
Standard Electrica, dimensionada 
ao âmbito mundial, ocupa uma 
prestigiosa posição nos mercados 
estrangeiros. À eficiência de uma 
administração actualizada, a capa- 
cidade técnica dos seus engenhe:- 
ros e o profissionalismo com que 
todos os empregados encaram as 


suas responsabilidades, são valo- 
res reais e positivos que, de ma- 
neira decisiva, contribuem para 
o engrandecimento da Standard 
Eléctrica e da indústria nacional 
Associada da International Tele- 
phone and Telegraph Corporation, 
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de 200 fábricas localizadas em 
mais de 60 paises, a Standard 
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C. D. U. 517.96 


THEORETICAL FOUNDATIONS OF THE FINITE 
ELEMENT METHOD 


RESUMO 


O metodo dos elementos Anitos é hoje o instrumento maus 
geral, e um dos mais eficientes, da análise de estruturas. 

E por outro lado uma técnica matemática geral, e o fim 
principal deste trabalho será apresentado como tal, utili- 
cando a Análise Funcional como linguagem ideal para a 
formulação abstracta geral, 

E fundamental para quem aplique o método saber prever 
a convergência para a solução exacta de uma sucessão de 
soluções aproximadas obtidas empregando elementos cada 
vez mais pequenos, 

No caso de se conseguir conformidade entre elementos, 
o método dos elementos finitos é um caso particular do de 
Ritz, de modo que pode assegurar-se a convergência desde 
que se corrija a complexidade. 

Um eritério geral de complexidade é pois justificado 
neste trabalho. Tal eritério exige que as componentes do 
campo, e todas as suas derivadas de ordem não maior que 
a da derivada de maior ordem que intervenha na expressão 
da densidade da energia de deformação, possam tomar 
qualquer valor constante dentro do elemento, 

Prova-se finalmente que um tal critério é também um 
critério de convergência, isto é, uma condição suficiente de 
convergência, mesmo em casos em que a conformidade entre 
elementos não se verifique. 


NOTATION 


All thé symbols will be defined where they are 
introduced. 

The following general iconvention are adop- 
ted: 
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SYNOPSIS 


The finite element method is nowadays the most general 
and one of the most powerful tools for the analysis of 
structures, 

dt is also a general mathematical technique and the 
main concerne oj the paper is to present it in this light. 
Functional Analysis is used as the ideal frame for a gene- 
ral abstract formulation, 

The ability to predict convergence to the exact solution 
of a sequence of approximate solutions obtained from 
patterns uf finite elements with decreasing size is funda- 
mental in the application of the method, 

In case conformity between elements is obtained, the finite 
element method is a particular case of Ritz'a method, so that 
convergence can be ensured as tar as completeness is achieved, 

À general completeness criterion is justified in the paper. 
Such eriterion requires that the field components and all 
their derivatives, of order not higher than the highest order 
of derivative entering into the energy density expression, 
can take up any constant value within the element, 

It is finally proved that such eriterion às also a general 
convergence eriterion, i.e. a suflicient conaition for conver- 
gence even if conformity is not achieved, 


1. Matrices (or vectors) will be denoted by 
bold face symbols: q, H. 
2. The dummy index convention will be used : 
= A,;%g + Agão T e 
3. A derivative will be denoted by a comma 
followed by indices indicating the variables with 


A 
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respect to which the function is differentiated. 
The order is indicated by a superscript in paren- 
theses : 


e em + 


4. À sequence will be denoted by its general 
term between braces: (U,). 


1. INTRODUCTION 


The finite element method is nowadays the 
most general and one of the most powerful tools 
for the analysis of structures. 

Although it was developed for structural anal- 
ysis it is really a general mathematical technique, 
and the main concern of this paper is to present 
it in this light. Mikhlin's book [1] was used as 
a basis for such purpose. 

The important problem of the convergence to 
the exact solution of a sequence of approximate 
solutions generated by elements with decreasing 
size will be given special attention. Experience 
seems indeed to indicate that the best control of 
the approximation error consists in examining 
the behaviour of a sequence of that kind. It has 
also been observed that no reasonable approxi- 
mate solutions are likely to be generated if the 
type of finite element used is such that conver- 
gence to the exact solution is not obtainable. 

Before convergence to the exact solution was 
given the attention it deserves, there was a ten- 
dency to make monotonic convergence play the 
fundamental role. Monotonic convergence is 
nowadays no more considered so important. It 
has indeed been demonstrated [2] that conformity 
(a condition for monotonic convergence) does 
not always speed up the convergence to the 
exact solution, i.e. less approximate solutions 
have been obtained for some problems when the 
monotonic convergence requirements were vere- 
fied than when they were not. 

The capacity for convergence to the exact 
solution of some kinds of elements has been 
already examined in th case where continuity is 
preserved [3]. The author himself presented a 
first proof [4] of the known criteria [2,5] which 
is also valid for cases where continuity is violated. 

If continuity is not violated, the finite element 
method becomes a particularization of the 
classical Ritz method. This connection with the 
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Ritz method has been observed very often but 
very seldom studied in detail and explored. 

It is very important to notice however that, if 
continuity is violated, the finite element method 
is not a simple application of the Ritz method. 
A section of this paper will be devoted to 
demonstrating that convergence to the exact 
solution is still possible even in those cases. 


2. STATEMENT OF THE PROBLEM 


Let 4 be a linear bounded operator defined 
for a dense linear subset M of a real Hilbert 
space H, Assume the operator À to be symmetric 
and positive definite [1]. 

Let (u, v) denote the scalar product of two 
elements of H. Let ||u || denote the norm of an 
element in H. 

This paper is concerned with the solution of 
the equation 


Ai== f (1) 


that is, in the determination of the element u 
which the operator A transforms into f; u and f 
are elements of H., 

Equality (1) is meaningful if element « belongs 
to M. It is possible however that no element 
of M can correspond to an arbitrary element f 
of H; this is what is meant by stating that 
equation (1) can have no solution in M. 

It can be shown [1] that, if equation (1) has 
a solution, this will be unique. It can also be 
demonstrated that the solution of equation (1) 
minimizes the funcional 


F (u) = (Au, u) — 2 (u, f) (2) 


and conversely, that the element which minimizes 
Fin M satisfies equation (1). 

If A is positive-bounded-below, and not merely 
positive definite, that is, if 


(Auu) =? (u,u) (3) 


y being a real constant, then the field of defini- 
tion of À can be extended so that equation (1) 
has a solution for an arbitrary element f of H. 

The extended field of definition belongs to a 
new Hilbert space, H ,, which is a dense subset 
of H, defined as the completion of the Hilbert 


TECNICA N.º 399-400 


space which results from associating with the 
elements of M the scalar product 


[4,0] = (Au,0). (4) 


This scalar product, which will be denoted by 
square brackets, is called the energy product. The 
norm in H, is termed the energy norm and will 
be denoted by bold vertical rules: 


Iul==V [u,u). (5) 
The energy norm of the difference of two 
elements is the distance between both: 
d (uv) = |ju—o|. (6) 
The square of the energy norm is termed 
energy [1]. 
If uy is the solution to equation (1), then 
Au, = fe (7) 
Functional (2) can thus be expressed as 


F(u) = [uu] — 2 [u,u]. (8) 
It can further be transformed into 


E (u) == [(u — ug), (u— ug)] — [ug , ug] = 
== uu —|u |. (9) 


Expression (9) makes it clear that the minimum 
value of Fin H, is obtained for u = u,. 

A sequence of elements (u,, |, belonging to 
the field of definition of a functional F, is termed 
minimizing [1] for F if 
) =H (10) 


um 


lim F(u 


Ho 


E, being the exact lower bound of F. 


As 
E = F(u)=— |u| (11) 
equation (10) implies 
lim | Han ÉS Ho | == 0. (12) 


Equation (12) means that any sequence which 
is minimizing for F converges in energy to the 
exact solution. Energy convergence is character- 
ized by the fact that the distance between each 
term of the sequence and its limit tends to zero /1.. 
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3. PARTICULARIZATION TO VECTOR 
FIELDS 


Let 2 be open, connected and bounded domain 
with a finite number of dimensions. Let S be its 
boundary, which is supposed to be composed by 
a finite number of closed, smooth or piecewise 
smooth stretches. 

Let 2 be the closed domain resulting from the 
combination of Q and 5. 

Take for space H the space of the real vector 
fields (with a fixed number of components) whose 
moduli are quadratically summable over 4. The 
scalar product of a pair of elements, « and v, will 
be given by the Lebesgue integral: 


(um) = f vao + f u'vds= [ u'vao 
t) 5 q 


(13) 


u and v being column vectores containing the 
components of « and v. 

The number of components of the vectors is 
independent of the number of dimensions of the 
domain. 

Equation (1) can be written in a more explicit 
form: 

Au = f (14) 


A being a matrix of operators. These operators 
are from now on assumed to be differencial. 

The fields beloging to M are not supposed to 
satisfy all the boundary conditions of the problem. 
Those which are necessarily satisfied by every 
field in M and by each field in H,, are termed 
principal boundary conditions. The remaining ones 
are called natural boundary conditions. 

Any field belonging to M is supposed to meet 
homogeneous principal boundary conditions. 
Besides, both the field and the derivatives invol- 
ved in 4 must be continuous. These derivatives 
will not however generally be continuous for 
every field in H,. 

The energy product between elements belon- 
ging to M can be computed by the use of equa- 
tions (13) and (4): 


[1,0] = [. (Au) v dº (15) 
Q 


An energy product involving elements in H, 
not belonging to M can be computed as the 
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limit of the energy product of a sequence of 
pairs of elements belonging to M. 

It is assumed that the expression (15) for the 
energy product can be transformed, by suitable 
partial integration, into 


[E em É (Ru)! L (Rv) dO, (16) 


|. is a square, symmetric and definite positive 
matrix, R a diferential operator. 
The energy of any element vin H , is given by 


fai] sm Í, (Ru L (Ru) do. (17) 


The expression under the integral sign receives 
the name of energy density. 

Assume that Ru involves derivatives of compo- 
nent w; with order not greater than p,. The de- 
rivatives of order (p,— 1) or less are termed 
principal derivatives. 

As the energy [u,u] of any field u belonging 
to H, must be finit, (Ru) has to be bounded 
almost everywhere in £!! for every field in H,. 
The field components and their principal deriva- 
tives must thus be continuous almost everywhere 
in O, 

In what follows, f will be supposed such that 
the exact solution falls into the subset C, CH, 
of the fields whose components are continuous 
everywhere in !), together with their principal 
derivatives. These continuity properties will be 
referred to in the text as principal continuity condi- 
tons. 


4. APPLICATION TO LINEAR THEORY OF 
ELASTICITY 


Elastic theories involve three kinds of magni- 
tudes: and whose 
vectors will be denoted by 7, : and u. 

These magnitudes are related by three kinds 
of field equations which can be symbolized as 
follows : 


stresses, straims displacements, 


(a) Equilibrium equations : 


Eo=X (18) 
(b) Strain-displacement relations : 
Du=: (19) 
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(c) Stress-strain relations: 


o = He. (20) 
E and D are differential operators, X is the 
vector of the body force density components, H 
is a symmetric positive definitive matrix. 
Equations (18), (19) and (20) are valid on L. 
On the boundary 5, the equilibrium equations 
become : 


No =p. (21) 

N is a matrix whose elements depend on the 
orientation of the normal vector at a given 
boundary point. p is the vector of the tractions 
applied to the boundary. 

The analysis of the equilibrium of elastic bodies 
reduces to finding the solution of the system of 
field equations (18), (19) and (20) whiçh satisfies 
certain boundary conditions. The simplest and 
most important types of boundary conditions can 
be expressed directly in terms of displacements 
or tractions applied to the boundary. Let S, and 
S, denote the portions of the boundary where 
tractions or displacements are respectively pre- 
scribed. 

Operators E and D and matrix N are such that 
the following relation holds if u is continuous: 


“cT(0u) dO= [ (ETud OA [ (N9)” 
PRMOE [E Tua 2 + | (No)/uas 
(22) 


In this relation, vectors 7 and U are not neces- 

sarily related by the stress-strain relations (20). 

Combining equations (18), (19) and (20), we 
obtain : 

EHDu =X. (23) 

Combination of equations (19), (20) and (21) 
vields : 

NHDu = p, (24) 


Comparing equations (23) and (24) with equa- 
tion (14), there results: 


cdi for points in! (25) 

f=X ' E 

= di | for points on 5. (26) 
= 
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It can easily be shown that operator A defined 
by equations (25) and (26) has the properties 
which were indicated in Section 2, if the dis- 
placement boundary conditions are enough to 
eliminate rigid body motion. 

Using equation (15) we obtain: 


[4,0] = [. (EHDu) vdO + f. (NHDOu)'  vds. 


(27) 


Equation (22) allows the transformation of (27) 
into: 


[u,v]= f, (Du) H(Dv) dl, (28) 


Operator R coincides thus with D and matrix 
L with H.* 


Functional F becomes: 


=| 


a Sl 


TH eao 2 f HudO > f pTuds 
t ; 


5 


(29) 


i.e. twice the total potential energy, if the dis- 
placement boundary conditions are supposed to 
be homogeneous, 

The theorem of the minimum total potential 
energy, which states that the exact solution is 
the one, from all the compatible elastic fields, 
which makes the total potential energy a mini- 
mum, is thus a particularization of the theorem 
which affirms that the solution of equation (1) 
makes F a minimum in the space of the fields 
with finite energy. 

The formulation which has been presented is 
quite general as it is valid not only for linear 
two and three-dimensional elasticity but also for 
linear theories of plates, shells and beams. 


* The equation 


K (Ru) TL(Ry) dt = fa (Au)! v do + IJ. (Au)T v ds 


which results from the combination of (15) and (16), per- 
forms thus in the general theory of the role of the work 
equation (22). 
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In the case of a plate, for instance, vector U 
contains the transverse displacement and two 
rotations, vector E contains the curvatures and 
the transverse shear strains, vector 7 contains 
the bending and twisting moments and the trans- 
verse shearing forces. 

Operator D involves derivatives of the first 
order. The principal derivatives are thus of order 
zero. This means that the elements of C, are 
elastic fields with displacement components con- 
tinuous everywhere in O. 

The principal boudary conditions, which are 
supposed to be homogeneous, are those invol- 
ving linear combinations of the displacement 
components. The natural boundary conditions are 
expressed in terms of stresses. 

A very frequent simplification in the analysis 
of plates, shells and beams consists in neglecting 
the transverse shear deformation. 

This makes it possible to reduce the number 
of the unknowns to one (the normal displace- 
ment) in the theory of plates, and to three (the 
normal displacement and the tangential displa- 
cements) in the theory of thin shells. 

The simplified theories represent by themselves 
also a particularization of the general problem 
formulated in Section 2. The only field components 
are now the independent unknowns [1]. 

The rotations become in the simplified theory 
first derivatives of the normal displacement. The 
corresponding energy density involves thus first 
derivatives of the tangential displacements and 
second derivatives of the normal displacement, 
so that the principal derivatives are the derivatives 
of first order of the normal displacement and 
the derivatives of order zero of the tangential 
displacements. 

The energy will be finit if the normal dis- 
placement and its first derivatives, as well as 
the tangential displacements, are continuous 
everywhere in 2. As the first derivatives of the 
normal displacement are the rotations, the elements 
of C, are still the elastic fields with all the 
displacement components continuous everywhere 
in O. 

Similar conclusions could easily be derived 
for beams. 

The principal continuity conditions are thus 
the same both in the simplified theories and in 
the corresponding theories where the transverse 
shear deformation is not neglected. 
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5. EQUIVALENT PROBLEM 


Consider the domain “! subdivided into a 
number of subdomains, 2!, O 0º ....,9,... 
and let H' be a real functional Hilbert space 
whose elements are fields defined on the general 
closed subdomain “”. The scalar product between 
any pair of elements, u” and v”, belonging to H' 


is defined by 


€ e e! a E o! e O ai 
(0) = [ ú vao + | nv ds = 


à o 


= u v dl. (30) 


Let H, be another Hilbert space (index n refers 
to a certain degree of subdivision of “ into 
subdomains). Each element u, GH, may be 
regarded as a piecewise defined field. It represents, 
however, not truly a single field defined on O, 
but a set of fields u (one per subdomain), 
belonging to the different spaces H”. Such fields 
are called subfields of u,. 

The scalar product in H, is defined by 


od = PD (0) (31) 


where u” and v” are the subfields of u, and »,, 
and 2 denotes a summation over the whole set 


of subdomains. 

Let M' be a dense linear subset of H”. Every 
field in M” is assumed to coincide, on the closed 
subdomain “4”, with an arbitrary field of M. 

Let M” be the field of definition of a linear, 
bounded and symmetric differential operator, 
defined to be such that, within £ 


Au = Au (32) 


and, on 5, 


| (Au) v dO= | (Ru) L(Rv) dO — 
- sf , nº 
— À (Au) v dO (33) 
“ Ta 


u ond v being arbitrary fields belonging to M”. 
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Make 


[1uº,o"] = [ (Ru')'L (Av) dO = 


(1 


e (Au) v dO = (A'u,o'). (34) 


Call M, the dense linear subset of H, whose 
elements have subfields belonging to the sets 
M”, and consider a linear subset of M, such 
that, given any pair of its elements, u, and v,, 


with subfields w and v, 


2 tio] 


can be properly chosen as their scalar product. 
Let H,. be the completion of the Hilbert space 
obtained by associating the scalar product (35) 
(termed energy product) to such subset. 

The distance between any two elements in 
H ,e will be given by 


the magnitude 


b,, Ed D,), = (35) 


= 8, 


d,, (mu, V,) mes V Lu, SE a (1, 


unas — 7) 
HM, Caiu 


(36) 


1 
whose elements are such that their subfields u”, 
corresponding to adjacent subdomains, take equal 
values at points lying on the common interface. 

Let H, and H, be the field of definition and 
the range of an operator B, such that, if 


Consider now the linear subspace of H 


"n” 


Bu =u, (37) 
the sum of the values taken on the interface of 
two adjacent subdomains by the corresponding 
subfields of u, is equal to the sum of the same 
values respecting «,. The effect of the operator 
B, is thus to distribute that sum equally between 
the subdomains in contact. 

Let v, belong to H, and let v” 
subfields. We can write 


denote its 


“ 


, Ea ER 1. 
(14, Dada == 2, (u”, U y cr 5 | o uv d 1 == 
e e 
à Í e T e ve 
=D f uivas= 3 urwry = 
e Es e 
Rs ( H, ! D, ), po ( B, H,, é b, ),, E (38) 
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Let M, be the field of definition of a new 


operator. 4,, such that 

A, H,, cs B, A (39) 
in which f, is an element of H, with subfields 
ff = Au. 


By virtue of equations (38) and (39), we have 


(A, Ugo 0) ani (B, fui ula ni (fis Di), era 
ge ba (AC u,0) == b» (MT) == [uns Value (40) 


Consider now the subset of H, whose elements 
fulfill the condition that all the corresponding 
subfields coincide, within their respective subdo- 
mains, with a given field u belonging to H, (the 
same for all the subdomains). It is clear that 
such subset is contained in H,. Call H, the 


Hilbert space obtained by associating the scalar 
product [1,,v,], to such subset As 


Lu, 2), há b> [w, DJ res Ef uy » ar dt 
(41) 


= | (Ru) L(Rv) dO = [4,2] 


[u,. v,]), is a proper scalar product. No contra- 
diction is thus introduced if H 4, is assumed to 
be contained in H ,.. 
Cali C,, the subset of Hy, corresponding to 
the subset C, of H . 
Consider the equation 
Au, ==, 


Ho mn 


(42) 


The solution of equation (42) in M, is gener- 
ally not unique and the operator 4, not positive 
definite. If, however, the field of definition of 4, 
is restricted to Ha,, then the operator A, 
becomes positive definite. Indeed, as Ha, is 
contained in H,, (40) permits to write 


(A, 4, u,, a ici Lu, H,, J, > 0. (43) 
Consider the functional 
> (1u,,) to [u, , a R dã 2 (Si EL (44) 
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Let u,, denote the solution of (42) in Ha, Then 


A n Hon f, ' (45) 
If u, belongs to H/ then 
(14, fo) pes UM A, Moss da e [u,, Ugo ne (46) 


Introducing (46) in (44), there results 


Edu) = uu, —2 [U,Ug = (4, — Ugo), 
(u,— u5,)] — [ge Ugu de (47) 

Expression (47) makes it clear that u,, mini- 
mizes F, and that any element which minimizes 
FE, in H, must coincide with u,,. The solution 
of (42) in H, is thus unique. 

Assume now that f, vanishes on the subdomain 
interfaces, and let f be any field in H which 
takes the same values, within the subdomains O”, 
as the subfields of f,. Let u, be the element 
in H,. which corresponds to the element u of 


H, Then 


ld =Duiry=D [ ut trad= 
e cv mn 


(48) 


-5/ tao + u'fds = (u, f). 
" tie Ss 


By virtue of (41) and (48), 


E (u,) = [u,u])—2 (u, f) = F(u). (49) 
Let u, be the of H, which corresponds to u,,. 
By virtue of (9), 


E (ug) = [ug — uo), (ug — u9)] — [ug 49]. (50) 


As F,(u,)F(u) and u,, minimizes F (u,) in 
Has uy minimizes F(u) in H,. But equation (50) 
shows that F can only be minimized by u 
if w, coincides with u,. Thus the solution of 
equation (42) in H, coincides, within each 
subdomain, with the solution of equation (1) 
in H,. 

This means that the problem of the solution 
of equation (42) is equivalent to the problem of 
the solution of equation (1). 
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The concepts introduced along this Section 
represent a generalization of the concepts intro- 
duced in the preceding ones. It is convenient to 
interpret such generalization in terms of three- 
dimensional Elasticity. 

Operator A has the same meaning for 
subdomain ! as operator 4 for the global 
domain, “. While A is defined by expressions 
(25) and (26), respectively for points located 
within 2 and on 5, operator A” is defined by 
the same expressions for points located within 
9º and on S. 

f. represents an external force distribution 
acting on each subdomain of the body. Assu- 
ming that f, belongs to H, is equivalent to 
assume that the total force acting on the interface 
between two adjacent subdomains is equally 
distributed between both. 

The value taken by f, on the subdomain 
boundaries are thus half the surface density 
values of the forces acting on those interfaces. 
These values must vanish, if the body force 
volume density is to be bounded everywhere 
in 2, 

To solve equation (42) means to determine an 
elastic field which verifies equation (23) within 
each subdomain and equation (24) on 5, and 
whose stresses present, on the subdomain 
interfaces, the discontinuities required to equili- 
brate the external forces applied on such inter- 
faces. Any solution of equation (42) equilibrates 
thus the external force distribution symbolized 
by f: 

The solution of equation (42) becomes also 
compatible if it belongs to Ha,, because the 
displacement boundary conditions are then 
respected on S, and the continuity of the 
displacements is preserved across the subdomain 
boundaries. The corresponding subfields coincide 


thus, within each subdomain, with the solution 
of equation (1). 


6. THE FINITE ELEMENT METHOD 


The finite element method is a general tech- 
nique of numerical analysis which provides an 
approximate solution for equation (1). 

In this method, domain “ is considered to be 
decomposed into a finite number of subdomains 
and families of fields are considered which have 
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different analytical expressions inside each 
subdomain. 

A finite element is a closed subdomain, 2”, 
together with the family of fields which are allo- 
wed to occur within it. This family is a linear 
combination with coefficients q; of a finite 
number of unit modes, so that each field of the 
family corresponds to ascribing particular values 
to the parameters q. 


The values of the field components and its 
principal derivatives, at a certain number of 
points on the boundary of the elements, called 
nodes or nodal points, are as a rule chosen as 
parameters. 

The type of an element refers to its general 
shape, nodal point specification and to the 
allowed fields, analytically defined by expressing 
a general field w in terms of the parameters 
and the coordinates with respect to a given frame: 
and the coordinates with respect to a given 
frame : 


= (rox, «JW (51) 


q is the vector of the parameters. * 

Elements P/, of matrix P* are supposed to be 
continuous and have continuous derivatives of 
order (p;— 1), or less, in the closed domain Qe 
occupied by the finite element e. The unit modes 
are definede by the columns of q“ 

We suppose that each field component depends 
only on its own values at the nodes and on the 
values taken by its principal derivatives also at 
the nodes. Thus, if q” corresponds to the field 
component w or one of its derivatives at any 
node of the element, all the magnitudes 21 ; for 
which k &i will be equal to zero. 

It q; corresponds to a derivative of order s of 
ui, 97; Wil take te form 


, | “m Xy N 
5, (xy, Xp» ..4 1º) ER (1 YS, e) 


* The allowed fields need not be introduced by 
giving the expression of the field components directly 
in terms of their own nodal values and the nodal 
values of their principal derivatives. They can indeed 
be given in terms of equal number of arbitrary para- 
meters which in turn can be expressed in terms of 
those nodal values (see |9]), 
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in which |* is a typical dimension of the element, 
for instance its maximum diameter, and Yí. is a 
function which does not depend on the absolute 
dimensions of the element. This is necessary in 
order that equation (51) can be homogeneous. 

The different finite elements are compatibili.. 
zed through the specification of reduced continuity 
conditions. These require that the values of the 
field components and their principal derivatives 
be the same at coincident nodes of adjacent ele- 
ments and equal the prescribed ones at the nodes 
located on 5S,, the portion of S where the prin- 
cipal boundary conditions are specified. 

A point of the domain is said to be a node of 
the system if it is a nod for one or more elements. 

Let q, be the vector of the field components 
and their principal derivatives at every node of 
the system but those which are located on 5. 
The reduced continuity condutions can be expres- 
sed by writing for each element 

qG=T], (53) 
where matrix T' depends on the topology of 
the system. 

Equations (53) show that the knowledge of q, 
is enough for the definition of the field within 
every element of the system. 

The reduced continuity conditions are gener- 
ally not sufficient to make the field components 
and their principal derivatives continuous across 
the element boundaries. This depends on the 
type of the element. 

If the type is such that the reduced conti- 
nuity conditions are sufficient to ensure conti- 
nuity of the components and their principal 
derivatives across the element boundaries, the 
piecewise defined fields generated by the system 
of finit elements are said to be conforming. Every 
conforming field thus belongs to C,,, i.e. to the 
subset of H, corresponding to C,. 

It the continuity requirements are violated 
across the element boundaries, the fields are said 
to be non-conforming. 

Let Ll, be the subset of H, containing the 
elements whose subfields are defined by equation 
(51) and compatibilized through the reduced 
continuity conditions. Any element u, Gl, can 
be expressed, withiin Q”, by 


=P" q, (54) 
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where 
O =e" 7º. (55) 
Take for space Ha (see Section 5) the space 
spanned by Ll, and Hy,. The distance between 
any pair of elements belonging to H ,. is defined 
by expression (36). The discussion of complete- 
ness and convergence will be based on that 
concept of distance. The distance between any 
element u, in LU, and any element «in H, will 
indeed be measured by the distance between u, 
and the element in H, corresponding to u. 


The approximate solution, u,,, which the finite 
element method -provides for equation (42), and 
thus for equation (1), is determined by making 
the functional F, stationary in U, . Such solution 
could be the exact one if up, was contained in 
U,. As, generally, it is not, the solution yielded 
by the finite element method is only approximate. 

Introducing (51) and (35) in (44), we obtain: 


E = Sgikq 2" 0") (56) 
where 
Ke = ) (Ro')L (R$ dO (57) 
tu 
= [ «tao. (58) 


4 
Introducing now (53) in (56), we obtain 


= (S TOK TI, — 29, (STO 05. (59) 


Making 
K= ST KT (60) 
O, = 51" 0 (61) 
there results 
- T T 
Fa ci q, K, q, o. 24, 0, . (62) 


The stationary conditions for F, are obtained 
by equating to zero the derivatives of F, with 
respect to the mutually independent parameters 
q,;. It results in the system of linear equations 


(63) 
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Introducing (57) in (60) and using (55), we 
obtain 


K => ) , (ROL(RO) dO (64) 


0, = 2, f pt dO, 
d Ta 


Matrix K 


»”M 


(65) 


is non-singular whenever the 


columns of P” are linearly independent. As L is 
definite positive, K, is also definite positive. 

If K, is non-singular, the parameters q,, can 
be uniquely determined by solving the system 
of equations (63). Let q,, be the vector of the 
parameters which verify equation (63). Functio- 
nal F, cam be expressed as 


F=0, K,9,—20, Ki, =(4,— 0,0)7K, 
(q, — To) cê Ko, (66) 


As K, is definite positive, the first term in 
the right-hand side of (66) is positive unless 
q, equals q,,. This proves that the solution of 
(63) minimizes F, in U,. 

Let now «,, and u,, be two elements belonging 
to U,. Let q,, and q,,, be the vectors of the 
corresponding parameters. The energy product 
of u,, and u,, can be given by 


Lu Ju" Uau) 4 == 3 f (Ru,) “L(Ruj) dO = 
e JD 
= 42 | (RP) URI) AO, 
e De 
T T Tr «Tera «TE 
o q in RT an a q,.0, pai Th, 3 /, z kr dt 


À E; 
na DI (a) ty da= BD [ uitçao 
e (e eu + i 
= uti) = (5, foda fiada 
(67) 


where f,, and f,, are the right hand sides (of 
equation (42)) which u,, and u,, correspond to 
(as approximate solutions). 

It results from (67) that the functional F, may 
take in LU, the following expression 
Eu, )=[u,,u),—2u,,f),=[u,,u,],—2[u RI: 


Hp an 


= [(4,— Uau) (4, — Un) a — [Man Man) n (68) 


which makes it clear that u, minimizes F, in U, 
Such expression will be used in Section 10. 


7. THE RITZ METHOD 


The method just described is justified if it can 
generate a sequence of fields converging to u, 
(the solution of equation (1)), when sucessive 
subdivisions are considered with elements of 
invariant type but decreasing size. 

Conditions to be met by matrix ? in order 
that this convergence may be ensured can be 
established if it is remarked that the finite ele- 
ment method is related to the well-known Ritz 
method (1, 6]. 

The Ritz method is a technique for generating 
a minimizing sequence for a given functional, 
say F. This technique. which can be used 
whenever H is a separable space, is based [6] on 
the determination of a sequence of families, (Va), 
satisfying the following conditions : 

(a) the sequence is complete in energy with 
respect to a class CCH , containing u, (comple- 
teness requirement) ; 

(b) the nth family depends on a finite number, 
N, of arbitrary parameters ; 

(c) every element which can be obtained by 
ascribing arbitrary values to the parameters 
belongs to C, (conformity requirement). 

In what concerns condition (a), it is remem- 
bered that a sequence of families of elements is 
said to be complete in energy with respect to a given 
class CcH, if it is possible, for a specified 
e >0, to find an integer N such that, in each 
family with order n > N, there exists an element 
u.« Which satisfies the inequality 


d(u,u,,) << E (69) 
where u is any element of €. 

The terms of the minimizing sequence ( Vu) 
are obtained by minimizing F in each family V, 

The elements of the nth family are generally 
given as a linear combination with coefficients 
qui Of N linearly independent fixed elements 'V,. 
which are termed coordinate elements : 


(70) 
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q, being the vector of the coefficients q,, and 'V, 
the matrix with columns 'F,. 

Family V, becomes thus a linear N-dimensio- 
nal space. Introducing (70) in (2), we obtain: 


F(u))=q,K,q,— 20,4, (71) 
in which 
se Í, (RETRO (72) 
and vector Q, is defined by 
Q = f ” 7 f dO (73) 


The values of the parameters which make F 
stationary can be determined by solving the 
system of linear equations 


K q, =0 (74) 


H 


K, is a non-singular matrix if the coordinate 
elements are linearly independent /1). The system 
has thus a unique solution which provides the 
unique stationary point of Fin V.. 

The finite element method can be considered 
as a technique for the application of the Ritz 
method only if the piecewise defined fields are 
conforming (conformity requirement). Only thus 
can indeed condition (c) be respected. The sets 
V, are then the subsets of H, corresponding to 
the subsets U, CH,. 

In order that convergence to the exact solu- 
tion may be obtained, it is thus only necessary 
to meet condition (a), that is completeness. It 
will be seen later on how this can be obtained. 

Comparing (70) with (54) it can be concluded 
that the coordinate fields used in the finite ele- 
ment method are defined by 

Wo=0"T'O=% within element e. (75) 

The analytical expression of the coordinate 
fields varies thus from element to element and 
this piecewise definition is the main characteristic 
of the finite element method. 

It is also important to notice that such piece- 
wise definition and the reduced continuity con- 
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ditions allow matrix K, and vector Q, to be 


assembled from simpler matrices, K' and 0”, con- 
nected with the finite elements themselves (see 
equations (60) and (61)). This is one of the 
most interesting features of the method. 


8. MONOTONIC CONVERGENCE 


Assume a sequence of families, (Va , fulfiling 
the conditions (b) and (c) stated in the preceding 
section and suppose that the nth family contains 
all the families with smaller order. As v,, makes 
F a minimum in V,, we have: 

E(D) > Fl) > EM) > ce > F (0) > cc E (Mg). 
(76) 


By Bolzano's theorem [7], the sequence [Fv,,)| 
converges to a limit which cannot be smaller than 
F(u,). It is remarked that this conclusion is valid 
even if condition (a) of Section 7 is not obeyed. 
If it is obeyed, then we know that the limit is 
E (ug). 

As the inequality 


F(v,)— F(v,,)<O (77) 
holds, for m < n, equation (9) yelds 
| Don — Ho | < | Dam — Ho . (78) 


This means that the distance to the exact 
solution decreases when n increases. Converg>nce 
is said to be monotonis. 

Monotonic convergence does not ensure con- 
vergence to the exact solution. On the other hand, 
convergence to the exact solution is not necessa- 
rily monotonic. 

Consider now a sequence of approximate solu- 
tions generated by finite elements with decreasing 
size. 

Conformity and the requirement that the 
family of fields corresponding to a given subdi- 
vision contains the families corresponding to 
elements with larger sizes have been proposed 
Melosh [8] as sufficient conditions for mono- 
tonic convergence of such sequence. 

However, as the approximate solution minimi- 
zes F. in LU, regardless of conformity being 
respected, the requirement that each family of 
fields contains the families corresponding to ele- 
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ments with larger sizes stands by itself as a 
sufficient condition for convergence. We can 
indeed write the set of inequalities (76) once this 
conditions is fulfilled 


9. COMPLETENESS CRITERION 


In what concerns convergence to the exact 
solution, we know that the Ritz method generates 
a minimizing sequence and that a minimizing 
sequence actually converges in energy to the 
exact solution. 

Convergence to the exact solution can thus be 
ensured if conformity and completeness are both 
achieved. We shall see however that complete- 
ness is the truly important requirement. 

Before proceeding further we remark that 
completeness of a sequence of families with 
respect to a set C CH, has a meaning provided 
we can computte the distance between every field 
of each family and any element in C (see Sec- 
tion 6). 

A general criterion for completeness will be 
stated and justified in this section. This criterion 
was presented in a recent book /9| by Zienkiewicz 
but is has not yet been justified as far as we 
know. 

Let (p, — 1) be the maximum order of the prin- 
cipal derivatives for component u.. 

We wish to demonstrate that completeness 
will be obtained if the general analytical expres- 
sion for n;, within element e (ses equation (51)), 
is given* as a polynomial with a number of 
arbitrary coefficients equal to the number of unit 
modes corresponding to the element. Furthermore 
this polynomial expression must contain a com- 
plete polynomial of the p,th degree all the terms 
of which are affected by independent arbitrary 
coefficients. The terms of higher degree can 
vanish whatever the values taken by those coef- 
ficients. 

We remark that, if this is the case, the field 
component «, or any of its derivatives of order 
p; or less can take any arbitrary constant value 
throughout the element if suitable values are 
ascribed to the parameters. In order that the deriv- 
ative u Nm assumes an arbitrary constant value 


V in Q, it is then indeed only necessary that 


(*) See footnote at page 442. 
pag 
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the coefficient which multiplies the monomial 
(x, %» “Jinu, be equalto V/r(r— 1)...s(s—1)..., 
all the remaining coefficients being equal to zero. 

The right hand side of equation (1) has been 
constrained in Section 3 to be such that solution 
uy belongs to C,, so that the derivatives of order 
p; Of solution «, are bounded but not necessarily 
continuous. 

In the next Sections we assume furthermore 
that the exact solution falls into a subset of C,, 
C,, such that the derivatives of order (p, + 1) 
of the field component «, are continuous within 
each element. Discontinuities of the p.th and 
(p; + 1) th derivatives are still allowed at points 
which always remain on element boundaries as 
the size of the elements is progressively reduced. 

Let C,, be the subset of H, corresponding to 
the subset C, or H,. 

Any field u, belonging to C,, can thus be 


represented inside O” by the following Taylor's 
expansion of its subfield components : 


u; = u; (0) + u;(0) (x, — x;) +... + 
1º ati 
+ ST ipod (O). (> 15) ee (ro) + 
1 eis 
(DD! Cite 
0 


«(x,—x,) 


+ 


HO) (x,—2;) (1,20) +.« 


(79) 


O and O, are points of $. O, depends on the 
coordinates of the point where u. is to be deter- 
mined. 

Let us consider now a polynomial field with 
components 


u; = u; (O) + u, ; (O) dx, — x;) |... + 
1 


+ pl o Sb (0) (x, — 5) 00 Kp— “) (80) 


which we call tangent field to « at O, 
As all the derivatives of order (p, + 1) are 


bounded inside !”, (79) and (80) yield: 


us ae 
to, + D! 


Va 


| º te | 
[M, ame 


(81) 


in which V, is an upper bound for all the (p, 4-1) th 
derivatives and [ is the maximum diameter of 
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element e. d is the total number of the (p, 41) th 
derivatives. 

By considering similar expansions for the de- 
rivatives of u,, it is possible to derive the following 
general inequality : 


etr) “a (1) d 


i-rt+l 
ad O hem Unid E GEE. (p,— r+D! Etr 


(82) 


for r< p,. 
As operator R involves derivatives of wu; of 
order p, or less, we have: 


[R(u — nO) LIR(u uv )<Vo(1)” (83) 


for  sufficiently small. V, is a positive number. 
Thus: 


E aó á N<a V, (190, (84) 


If tangent fields wu” are considered for every 


subdomain “”, piecewise defining a field min 22, 
we obtain, by using (35), 


(|u, Ni mc ao = [(u—u ), (uu, )L< Vol 


(85) 


in which |, denotes the maximum value of |" in 
the whole set of elements. 

This means that the distance between any field 
in C, and the tangent field u, , piecewise defined 
by (80), tends to zero with the size of the finite 
elements. 

Consider now a type of finite element gener- 
ating a sequence of families of fields whose 
completeness is to be investigated. 

Call «/” the field within the finite element e 
such that the values of its components and their 
principal derivatives at the nodal points are res- 
pectively equal to the values of the components 
of the field u,GC,, and corresponding partial 
derivatives at the same points. 

Suppose the general criterion to be satisfied. u” 
can thus be one of the fields which can occur 
within the finite element. Let this field correspond 
to values q; of the parameters: 


te 


Tue À Tg (86) 
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On the other hand 


vt =" (87) 
From (86) an (87) we obtain 
[nf us 5)! (88) 
or, considering (52), 
uy |=S 0. | di (aj — ai] (89) 
: 


in which s is the order of the field derivative to 
which parameter q, corresponds. 


But 
4 1º 4 +g , 
gv ( pe Td ) (90) 
and 
| a 
(7) má O O a (91) 
IPA cul try” dtxlT) ccoaltzdr | 
For this reason, 
ui ug | =| Paula; —a5)] 
E say RE é E 
j 9 (x)... d(x,/1) 
(92) 


As the absolute dimensions of the element do 
not appear explicitly in the functions Vis these 
functions remain bounded as the size of the 
element decreases. 

The same happens to the derivatives 0" 


Ldtxih)os. 

tions 7; and their derivatives of order (p, 1) 

or less were in Section 5 supposed to be con- 

tinuous. Assume the moduli of all these magni- 

tudes remain below a positive number V.. 
Then 


d (x,/1'), for r < p, because the func- 


| telri fer) 
dE Ru O 


Ee E O Roo 
(93) 


On the other hand, as the components of u” 
and their principal derivatives take the same 
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values at the nodes as the corresponding magni- 
tudes in u, equation (82) permits us to write: 


flo d E E soa a 
+ | (p;—s+])! Vth) 


| te 
4; —4 
(94) 


when parameter q, corresponds to a derivative 
of order s. 

Equations (93; and (94) hold even if the p,th 
derivatives of wu” are discontinuous in O”. This 
is an important remark because sometimes |10] 
the element itself is considered subdivided into 
parts and the field admits different analytical 
expressions within each part. Our proof remains 
valid however even if the derivatives of order 
p, are not continuous across the internal boun- 
daries of the element. 

As the parameters cannot correspond to deriv- 
atives of order larger than (p, —1), s cannot be 
larger than (p,—l) and equation (94) yields 


F te —qf A V, (1º) pt =s+1. (95) 
| 2! 


Introducing equations (95) in (93) we obtain: 
te(r) fetr) V, gepi-s+1 e 
Pads Rakul | < Vs “a (1) dN (96) 


Nº“ being the total number of parameters corre- 
sponding to element e. 

This equation is still valid for r==0, ifthe 
derivatives of order zero are interpreted as the 
field components themselves. 

The similarity between equations (96) and (82) 
allows a jump straight to the inequality : 


( | u . ee : N )2 4 V, [ NR) (97) 


V, being a positive number and w/ denoting the 
piecewise defined field which coincides with u” 
within a general element e. 

Equation (97) means that the distance between 
u', and wu” tends to zero with | . Combining (97) 
and (85) we conclude that the distance between 
u, and u” tends also to zero when the size of 
the element decreases, so that, as u, is an arbitrary 
element of C,., the completeness proof is finally 
achieved. 


fui 
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10. CONVERGENCE DISCUSSION 


Consider any type of finite element which can 
generate a sequence (LI, | of families of generally 
non-conforming fields complete in energy with 
respect to C,. 

We wish to investigate if the sequence of 
approximate solutions (u,,! obtained by mini- 
mizing F, in each family U, converges in energy 
to the exact solution. 

We know already that completeness implies 
convergence to the exact solution if it is associated 
with conformity. It will be concluded in this 
Section that completeness with respect to C, is 
a sufficient condition for convergence, regardless 
of conformity being obtained. 

Let w,, be the field in U, which presents the 
same nodal values of the field components and 
corresponding principal derivatives as 1,, (the 
solution of equation (42) in H ,,). As completeness 
is ensured, it is possible to determine N such 
that, for n > N, 


d, (Mia , Un) < E (98) 
: being a positive and arbitrarily small number. 
As F, is continuous, we can find : such that 
E (u,,) = E (um) + € (99) 
* being also positive and arbitrarily small. 
As u,, belongs to Ul,, and u,, (the approximate 
solution yielded by the finite element method) 
minimizes F in U,, 


F, (1,0) aa F. (14,1) (100) 

and 
E ua) < Eai E (101) 
Let now f,. be an element of H, defined by 
ar imeti A, Hu x (102) 
As u,, is an approximate solution to equation 


(42), f., generally does not coincide with f.. 
Let u,, denote the solution of the equation 


= f, (103) 


uu 


Au. 


in H,.. Assume that u,, belongs to C,,. This 
assumption will later be discussed. 
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Let u,, be the field in LI, whose components 
and corresponding principal derivatives take the 
same nodal values as u,,. As u,, belongs to C,,, 
and the sequence (LU | is complete with respect 
to C,, it is possible to find N, such that, for 
n>N,, 


d, (1, 4 Upn) E. er (104) 

&” being a positive number, arbitrarily small. 
As A, is a continuous operator, it is also 

possible to determine :” such that (104) implies 


| A, Hoy — A, 4 py | < Ee (105) 
&“ being positive and arbitrarily small. 
Let 
fim E A, Ubn p (106) 
As u,, is a solution to equation (103), 
Ay Mia = font (107) 


The inequality (105) can thus be transformed 
into 


| Fam et Fm | A” 4 (108) 


By virtue of (67), as u,, and u,, both belong 
to: E; 


ld, (un 1 Un Di Es [ua id Upy)o (ua E “bn MJ, (109) 
e ((u,, cá Upa) (Fa — fonDu . 


As, by Cauchy's inequality, 


(ua — Hu ), (em EM fon D, E | “an ra pn | 


Hm 


| Fan Fon IA a e Van "Ee 4 bn 4 (110) 
we obtain 
da (Manos) CV (E ua — Mou) (111) 


Combining (104) with (111), there results 


d, (am ? Un) < d, (tan , Un) se d, (u,, , Un) < iv 
(112) 


where 


Sette” Sua 


Lan Uby | 
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As É is a continuous functional, it is then 


H" 
possible, given : >> 0, to determine «” such that 


F (ua) = Fun) e. (114) 
As u,, belongs to ear 
Fo (a) > Fuga) (115) 
Thus 
Es (14.4) + e E a (4) . (116) 


Combining (99), (100) and (116), there results 


EM) E EE Pt do 
and thus 
Fi) = Eta) + 8º (118) 
in which "<<: + e”. 
But, as u,, and u, both belong to U,, (68) 
permits to write 
Eu) — F, (a) id [un ie Uns (144 ds Us )), ii 
e ld, (an Un]? = a (119) 


Combining (119) and (98), we obtain finally 


d Ma Un) AM Mia) E do Um) CE VEM, 
(120) 

Expression (120) shows that u, converges 
to Uop 

It remains to prove that u,, belongs to C,,. 
This assumption was indeed used to obtain (104). 

Our reasoning will be based on a theorem 
which is known to be valid for Poisson's 
equation 

Au = f (121) 

Such theorem states that u has continuos 
second order derivates in domain & whenever f 
is a Hôlder continuous function in “ [12]. 

A corresponding theorem is lacking which 
refers to the general problem with which the 
present paper is concerned. We are thus not sure 
that the derivatives involved in the operator 4 are 
continuous whenever f is Hólder continuous. It 
seems however very reasonable to expect the 
theorem to be true, at least as far as linear elastic 
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theories are concerned. It will thus be admitted 
that, at least in case of Elasticity, the Hólder 
continuity of the body force density implies the 
continuity of the displacement derivativesinvolved 
in the operator. 

u.. denotes the solution to equation (103) in 


eu 


H, . This means, in terms of Elasticity, that u,, 


represents the compatible field which equili- 
brates the same external forces as u,.. Such 
forces are of two kinds: body forces distributed 
within each subdomain and forces distributed 
on the subdomain interfaces and on S.u,, will 
belong to C,,, i.e. its (p, + 1)th derivatives will be 
continuous within “Y, if the derivatives involved 
in the operator are continuous within !”, and 
thus if the body force density corresponding 
to 4, is Hólder continuous within 2. 

The problem now consists in proving that the 
body force density corresponding to «,, is Hôlder 
continuous within S, no matter how large is n. 
We shall not attempt to investigate the general 
conditions in which such a statement is true. 

Sometimes, however, the proof is trivial. This 
is namely the case if the type of the element is 
such that the body force density vanishes or is 
forced to a prescribed bounded and continuous 
polynomial variation within each element, no 
matter the values of the parameters. 

If such a common situation arises, there 
remains no doubt that the completeness criterion 
is also a convergence criterion, even if conformity 
is not achieved. But, if it is achieved, the finite 
element method becomes a particularization of 
the Ritz method, and completeness will ensure 
convergence in any case. 

Our reasoning can be adapted to cases [10] 
in which the elements are subdivided into parts 
and the allowed fields have different analytical 
expressions within each part. The body force 
density is then generally not continuous within 
the element taken as a whole. Convergence will 
however easily be proved if each part is treated 
as a separate element. 


11. CONCLUSIONS 


The finite element method has been presented 
as an analytical technique which can be applied 
to a very broad class of problems. 
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Functional Analysis provides the frame for 
an abstract formulation in which some genera- 
lized concepts, like energy and distance, with 
a physical or geometrical origin, play a funda- 
mental role. Particularly, the definition of distance 
between two fields, as the square root of the 
energy of their difference, is an extremely con- 
venient basis for the discussion of convergence. 

The description of the finite element method 
also required the introduction of the concept of 
principal derivatives. Continuity of the field 
components and their principal derivatives is 
necessary if the energy density is to be finite. 
In two and three-dimensional Elasticity, for 
instance, as the continuity of the displacement 
components is enough to ensure a finite energy, 
the principal derivatives are of order zero. 

The fact that the finite element method is 
based on the decomposition of the global domain 
into subdomains, made it convenient to trans- 
form the initial problem into an equivalent one. 
In Elasticity, for instance, the initial problem 
consists in the determination of an elastic field 
which verifies the field equations everywhere 
within the domain, while the transformed 
problem consists in the determination of a set 
of fields respecting the field equations within 
each subdomain and verifying compatibility and 
equilibrium conditions on the subdomain inter- 
faces. In both cases the conditions imposed on 
the external boundary must be fulfilled. 

In the finite element method, the principal 
continuity requirements are replaced by reduced 
continuity conditions which may imply the 
fulfillment of the principal continuity conditions 
everywhere in the domain. If they do, the 
piecewise defined fields are said to be conforming 
but non-conforming if it happens otherwise. 

This is an opportunity to remark that con- 
formity has been with difficulty obtained for plate 
and shell elements [10]. Such difficulty results 
from the fact that rotations are usually regarded 
as derivatives of the transverse displacements. 
This has not however to be so (see Section 4), 
and conformity can be easily obtained if the 
rotations are considered as true displacements. 

When conformity is achieved, the finite 
element method becomes a particularization of 
the Ritz method. Such particularization is charac- 
terized by the piecewise definition of the field, 
which, together with the reduced continuity 
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conditions, allows matrix K, and vector Q, (the 
stiffness matrix and the force vector in elastic 
problems) to be assembled from simpler matrices 
and vectors which refer to each finite element. 

The matrix analysis which is developed in this 
paper for the determination of K, and Q, is 
a generalization of the displacement method of 
Stuctural Analysis. The force method [4] could 
also be used. 

Completeness is a sufficient condition for 
convergence to the exact solution in the Ritz 
method, i.e. if conformity is achieved. 

A general completeness criterion is justified in 
Section 9. Itis proved that such criterion ensures 
completeness with respect to a set C, containing 
the fields whose derivatives of order up to (p, + 1) 
are continuous within the subdomains corre- 
sponding to the elements. The principal deriva- 
tives are of order p, — 1, so that the criterion does 
not ensure completeness with respect to the set 
of all the fields with finite energy. 

In two and three-dimensional Elasticity, the 
(p; + 1)th derivatives are second order deriva- 
tives. Their continuity implies the continuity of 
the body force distribution density. The continuity 
of the body force distribution density is not 
however a very strong restriction, as disconti- 
nuities of the first and second derivatives of the 
displacements are still admitted at points which 
remain on element boundaries when the size of 
the elements is progressively reduced. The solution 
of problems in which external forces are distributed 
on element interfaces is thus not excluded. 

Completeness with respect to C, is a sufficient 
condition for convergence whenever conformity 
is achieved. However, it was shown in Section 10, 
that completeness implies convergence in any 
case, i.e., even when conformity is not achieved, 
if the body force density remains continuous 
and bounded within each element as n tends to 
infinity. 

Along the whole paper the principal boundary 
conditions (which correspond to displacement 
boundary conditions in Elasticity) were supposed 
to be homogeneous. If they are not homogeneous, 
the finite element method can however still be 
applied. All that must be done is to make the 
values of the field components and their principal 
derivatives coincide with the prescribed values, 
at the nodes which are located on 5,. As the 
size of the elements tends to zero, we obtain 
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approximate solutions tending to a solution which 
obeys the field equation inside the domain (in 
case the convergence criterion has been respected) 
and the prescribed boundary conditions on the 
boundary. This is of course the exact solution. 

It remains to indicate that the formulation 
presented in this paper is not the only possible one. 

In the present paper, the nodal values of the 
field components and their principal derivatives 
are indeed chosen as parameters in terms of 
which the reduced continuity conditions are to 
be expressed. 

It is however possible, in elastic problems, to 
take as parameters the resultants and moments 
of the forces distributed on the element boundaries. 
The analysis starts then from reduced equilibrium 
conditions, which are directly expressed in terms 
of such parameters [11]. The interest of this 
second formulation is that it can generate 
equilibrated solutions while the first formulation 
leads to compatible ones (if conformity is 
achieved). 

This second formulation may be generalized, 
as well as the first, to cover the general problem 
with which the present paper is concerned. It 
will be shown in a next paper how this can be 
done and which criterion may be used to ensure 
convergence to the exact solution. 
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APPENDIX 


COMMENTS ON THE PAPER: THEO- 
RETICAL FOUNDATION OF THE FINITE 
ELEMENT METHOD % 


|x rHE summer of 1965 following some clues 
in the literature [1] and also some early indica- 
tions in the experiments on the Zienkiewicz 
plate bending triangle [2], certain definite 
conclusions were reached concerning the conver- 
gence of non-conforming plate elements. The 
numerical experiments and the theory that 
subsequently reinforced the conclusions are 
described in the addendum to Ref. |2]. The 
experiments are now described less succinctly, 
as their import appears to have generally 
misunderstood. 

If we postulate that a particular solution 
converges as the mesh is refined, giving smoothly 
varying stresses, we may presume also that the 
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finite element representation gives sensibly 
constant bending stress within each element 
and hence sensibly constant curvature. This is 
why, in engineering terms, displacement functions 
must be “complete'|3). However, early experi- 
ments suggested that the responsabilities were 
not entirely fulfilled unless conformity was also 
imposed. It was fortunate that four elements 
were simultaneously under test |2], one of which, 
the Zienkiewicz version, was non-conforming, 
the other three conforming. All four gave the 
states of constant curvature under suitable nodal 
action and all four were derived from a single 
shape function routine with options, by numerical 
integration. In two cases this was exact. The 
routine was checked thoroughly by various 
numerical tests. 

The first indication of the strange behaviour 
of the nonconforming triangle was when mesh 
(a) of Fig. 1 responded correctly to distributed 
bending moment applied at the free end 
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16. 1, The mesh pattern on which experiments demon- 
strated convergence or non-convergence. 


whereas (b) did not. The evidence was already 
irrefutable, for the conforming triangles all 
behaved correctly in both mesh (a) and (b). This 
evidence suggested that a mesh of type (c), 
composed of three sets of parallel lines, gives 
convergence, whereas mesh (d), composed of 
“union jack” figures, would not give convergence, 
however fine the mesh. This suggestion may be 
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easily and decisively checked, by imposing a state 
of constant curvature on the basic figures which, 
replicated in both directions, give (c) or (d), 
Thus, when the external nodes of (e) and (f) are 
given slopes and deflections corresponding to 
states of constant curvature, either the elements 
behave, or they misbehave. That is, the internal 
node either takes up its correct slopes and 
deflection corresponding to the required state of 
constant curvature, or it does not. With the 
non-conforming triangle, (e) behaved but (f) did 
not. (Despite this result, the element is useful 
and gives accurate results in practice.) 

These observations did not complete the 
research. The theory they suggested showed that, 
provided the element has certain symmetries 
possessed by every practical element, mesh (c) 
converges, regardless of the angles of the three 
sets of parallel lines. This property extends also 
to the general anisotropic case, and it even 
extends to the case where the integration 
formula is not accurate enough to give correct 
stiffness coefficients [4]. As a bonus, the theory 
predicts correctly [1] that a mesh of equal Melosh 
rectangles converges. 

These observations and theoretical predictions 
seem to be at variance with the conclusions of 
Ref. [2]. It is the purpose of the present note to 
clarify what has been proved. 

The important point is that the convergence 
proof presented in Ref. [2] depends on u,, 
belonging to C,,, as indicated in the text. 

The analogy with the case of Poisson's equation 
permits the conclusion that u, belongs to C,, im 
the case when operator À is a second order 
differencial operator and the components of f, 
(body force density components in the elastic 
case) corresponding to the successive approxi- 
mate solutions, remain continuous and bounded 
within each element as the size decreases indefi- 
nitely. This is the case in two- and three-dimen- 
sional elasticity, and hense also for plates, shells 
and beams if the transverse shear deformation has 
not been neglected. 

The situation changes however if the order 
of the highest derivative of u, contained 
in A exceeds p + 1, because the boundedness 
of the components of f, does not then ensure 
the boundedness of the (p, + 1)th derivatives of 
the field component «. The (p, + 1)th deriva- 
tives can, for instance, be constant within each 
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element but unbounded as the size decreases 
indefinitely, whereas the (p, + 2)th derivatives 
vanish and are thus bounded. 

Such is the elastic theories which result from 
neglecting the transverse shear deformation. The 
appropriate convergence conditions can be es- 
tablished by resorting to the convergence condi- 
tions of the corresponding theories in which the 
transverse shear deformation is not neglected. 

In the simplified theory of thin plates, for 
instance, À is a fourth order operator and (p, + 1) 
is equal to 3. The convergence conditions are, 
completeness, and the condition that the third 
derivatives of the transverse displacement, corre- 
sponding to the successive approximate solutions, 
remain bounded within each element as the size 
decreases (see Ref. [5]). 

The union jack mesh (d) does not give con- 
vergence because the third derivatives of the 
transverse displacement are unbounded as the 
size of the element decreases. Such derivatives 
remain bounded in mesh (c): convergence is thus 
obtained. If conformity were guaranteed, com- 
pleteness would be sufficient for convergence 
regardless of the mesh pattern. Completeness 
and the boundedness of the (p. + 1)th deriva- 
tives corresponding to the successive approxi- 
mate solutions appear to be sufficient conver- 
gence conditions in the general case. 

Such conditions can at least be shown to be 
sufficient if a type of element exists with the 
same general shape and nodes as the one which 
is to be investigated, but meeting both complete- 
ness and conformity conditions. Let Un be the 
field generated by such an element, when the 
nodes and nodal displacements are the same 
as those corresponding to «,.. The same argu- 
ment that was used to prove the completeness 
criterion in section 9 of Ref. [3|) may indeed be 
used again (if the (p, + 1)th derivatives of u,. 
are bounded) to show that the distance between 
the two fields “a and «,. cam be smaller than 
a given positive and arbitrarily small number. 
As u,, is a compatible field, the inequality (116) 
of Ref. [3] remains valid and the convergence 
proof continues as in Ref. [3]. 
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PROGRESSOS RECENTES 


RESUMO 


Nesta conferência fez-se uma revisão de alguns dos pro- 
gressos mais prometedores na indústria do aço que poderão 
vir a afectar o equilíbrio entre a procura e a capacidade 
de produção Analisaram-se, portanto, técnicas recentes 
de redução do ferro, aciaria e fabricação no que diz respeito 
aos seus problemas tecnológicos e viabilidade econômica, 
Foi dada ênfase ao papel que o controle desempenha no 
«management» e na tecnologia. Discutiu-se ainda, sumãd- 
riamente, a importância que estes novos progressos poderão 
vir a ter no futuro da indústria portuguesa do aço, 


1. INTRODUÇÃO 


O aço é um dos mais baratos materiais de 
construção em engenharia. A despeito dos pro- 
gressos no abastecimento de outros metais e na 
produção de plásticos, os mercados em que estes 
competirão com o aço estão provavelmente limi- 
tados aos campos de embalagem e estruturas 
leves, e mesmo nestes a indústria do aço pode 
competir criando novos produtos tais como aço 
revestido a plástico que reune as vantagens des- 
tes dois materiais. 

Ao longo dos séculos a produção de aço 
aumentou sistematicamente, àparte períodos de 
depressão económica mundial como a dos prin- 
cípios da década de 1930 (fig. 1). Uma extrapo- 
lação da evolução da produção bruta mundial de 
aço feita por Lord Layton (!), director comercial 
de British Steel Corporation, sugere uma pro- 
dução total de 800 milhões de toneladas de aço 
por ano para 1985. Esta data foi tida como apro- 
priada para uma previsão a longo prazo. 

O consumo de aço nos países desenvolvidos 
está grosseiramente relacionado com o rendi- 
mento «per capita» médio do país, 30-50 kg/$100 
U.S./ano. Os Estados Unidos consomem cerca 
de 650 kg/cabeça/ano e a média mundial é apro- 
ximadamente de 150 kg/cabeça. Podemos então 
esperar um eventual aumento de procura em 
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Professor no Departamento de Metalurgia 
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SYNOPSIS 


Some of the most promising developments in the steel 
industry which are likely to affect the balance between 
demand and production capacity have been reviewed, Recent 
techniques used in ironmaking, steelmaking and fabrication 
have therefore been analysed with respect to their technolo- 
gical problems and economic viability, Special emphasis 
was given to the part that control plays in management and 
in technology. The importance that these new developments 
may have in the future of the portuguese steel industry was 
briefly discussed, 


Portugal da ordem de 5000 000 ton/ano. Pode- 
mos comparar este número com o verificado em 
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Fig. 1 — Crescimento da produção mundial de aço 
de 1890 a 1968, 


* Conferência organizada pelo Laboratório de Física e Engenharia Nucleares (Grupo de Metalurgia) e proferida 
no anfiteatro do Laboratório Nacional de Engenharia Civil em 25 de Setembro de 1969. 
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1967, que foi inferior a 1000 000 de toneladas, 
das quais uma fracção importante foi importada. 
Será então de esperar um substancial aumento 
no consumo de aço à medida que a prosperidade 
do país aumenta e talvez não seja insensato pre- 
ver uma taxa de aumento de procura de 8º/, por 
ano, até atingir 3 000 000 toneladas/ano em 1985.* 
A fig. 2 dá uma ideia do que isto significaria se 


«rails», dormentes e «fishplates». « . . 35000ton 
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Fig. 2 — Distribuição de uma produção de 3 milhões 
de toneladas de aço em proporções semelhantes às do 
U.K. em 1967. 


a distribuição dos produtos fosse a mesma que 
se verifica no Reino Unido. Contudo deve-se 
acrescentar desde já que esta não será na reali- 
dade uma comparação muito aceitável na medida 
em que os moldes de consumo variam de país 
para país. A intenção é, principalmente, dar uma 
ordem de grandeza dos produtos requeridos. 
Nos últimos dez anos, a capacidade de produ- 
ção excedeu em 50 milhões de toneladas a pro- 
dução mundial total; no entanto como este 
número apenas representa 10"/o da referida capa- 
cidade, ele não é certamente suficiente para 
afectar a viabilidade económica de uma siderur- 
gia moderna. Contudo, esta tendência para a 
capacidade produtora exceder ligeiramente a pro- 
dução representa uma ameaça para unidades 
obsoletas e para organizações usando práticas 
antiquadas. Há necessidade de olhar continua- 
mente em frente, o que implica a adopção das 
mais modernas técnicas de administração para 


* A taxa internacional é, presentemente, 8 a 12º/, por 
ano. 
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assegurar a «performance» e perspectivas de 
cada empreendimento. Algumas das técnicas 
essenciais para o estabelecimento da viabilidade 
económica dos actuais e dos novos processos 
estão sintetizadas nos pontos seguintes: 


ii Planeamento 


i.1 — pesquisas extensivas ao mercado 


i.2 — avaliação da rentabilidade dos inves- 
timentos usando técnicas de «dis- 
counted cash flow» 


i.3 — custos marginais para expansão da 
produção 


1.4 — «Network analyses» 


ii. Controle 


ii.1 — custos standard de produção 


ji.2 —técnicas de controle de «stocks» 
para minimizar os custos de arma- 
zenagem» 


ii.3 — utilização de técnicas de programa- 


ção linear 


ii.4 — controle de qualidade para garantir 
e aceitabilidade do produto. 


Estas técnicas devem ser usadas em cada caso 
particular para se estabelecerem objectivos defi- 
nidos a todos os níveis de administração. Infeliz- 
mente há poucas conclusões que possam ser 
transpostas de um local para outro e os proble- 
mas devem em cada situação ser abordados de 
novo. 

A escolha de uma nova unidade e de novos 
processos está também sujeita aos habituais con- 
dicionalismos metalúrgicos devidos ao tipo de 
produto requerido e ao número de unidades 
necessárias para prover a uma indispensávél uni- 
formidade do produto (muitas vezes, mais de 
50º, da produção deverá depender apenas de 
uma unidade). Este facto pode ter um efeito 
estratégico determinante se estamos entregues a 
determinado tipo de empreendimento e deseja- 
mos mudar para outro. 
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Nesta ordem de ideias pensamos que a indús- 
tria do aço em Portugal deparará, nas próximas 
duas décadas, com uma série interessante de 
opções. A procura expectável para o aço em 
Portugal é pequena comparada com a capacidade 
prevista para várias grandes siderurgias na 
Europa e Japão. Deste modo, as unidades mais 
económicas a instalar serão provávelmente dife- 
rentes das que seriam escolhidas para países com 
maior procura de aço. Isto é ma maior parte 
devido ao facto de que, como Leckie e Morris (2) 
salientaram, acima de 1,5 milhões de toneladas/ 
/ano, a rentabilidade de uma empresa depende 
mais da percentagem da capacidade de produção 
utilizada do que das dimensões das instalações. 
É provâvelmente mais económica uma expansão 
baseada em certo número de unidades de tamanho 
médio do que num menor número de unidades 
de grandes dimensões, porque assim se poderá 
seguir mais de perto as tendências do mercado 
de vendas. 

Propomo-nos abordar primeiramente a fabri- 
cação de gusa, depois a aciaria e finalmente o 
fabrico de produtos intermédios e acabados. 


2. PRODUÇÃO DE GUSA 


Devem ser considerados dois caminhos: 


1—o tradicional alto forno, produzindo gusa 
com alto teor de carbono. 


2 —técnicas de redução directa produzindo 
material com 92-95º/o de ferro, do qual 
92 a 97º/, no estado metálico, ex.: por 
um forno rotativo SLRN ou pela técnica 
de leito fluidizado ESSO FIOR. 


2.1 — Altos fornos 


Estas unidades desenvolveram-se grandemente 
em dimensões e produção nos últimos 20 anos. 
O maior alto forno japonês, FUKYAMA 2, tem 
um diâmetro de cadinho de 11,2m e produz 
6000 toneladas de gusa por dia (3). A mesma 
companhia está presentemente a construir um 
alto forno com 11,8m de diâmetro para produ- 
zir 7000 toneladas por dia. De acordo com 
Leckie (4), há uma relação geral entre os custos 
de funcionamento, produção e dimensões, que 
favorece as grandes unidades desde que estas 
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sejam utilizadas perto da sua capacidade nomi- 
nal, * 

As dimensões, por si só, não são garantia de 
uma produção volumosa ou económica. Em par- 
ticular é necessário um bom controle da marcha 
do alto forno. A marcha dos altos fornos japo- 
neses tem fornecido resultados excelentes em 
virtude da atenção dada às características a 
seguir mencionadas. Estes factores conduzem a 
uma marcha mais rápida com menor consumo 
de coque. 


1 — Controle da carga 


O alto forno Fukuyama utiliza : 
a) 70º/, sinter (6,7-50,8 mm, < 5'/, 
de O a 5 mm) de basicidade controlada. 


* É notável que as mesmas considerações se podem 
aplicar aos custos de frete e usando mineraleiros de 
100 000 toneladas estes custos para o Japão descerem a 
1/3 dos que se verificaram em 1960. 

A fig. 3 mostra como os custos descem ccm o aumento 
de produção. 


Segundo Leche 


so 


iToneiada 


Estudos 


Fig. 3 — Custos operacionais do A.F. em função da pro- 

dução (segundo Leckie) (4). Representa-se o excesso 

de custos relativos aos verificados num A.F. com 

9,75 m de diâmetro, 1 750 kg/ton de carga e B.O,]. = go, 
em função da produção. 


A fig. 4, por outro lado, mostra claramente os incon- 
venientes de trabalhar abaixo da capacidade óptima. 

Embora ao longo desta exposição se tenham conver- 
tido os custos ingleses em escudos, isto não significa que 
se considere que os custos avaliados tenham algum signi- 
ficado absoluto no caso português Eles são apenas dados 
para ilustrar tendências que têm um significado geral. 
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b) 30º/, de minério calibrado (9,5-25,4 mm) 
c) 475 kg de coque/tonelada, de alta resis- 
tência mecânica para manter a permeabi- 
lidade, com baixo conteúdo em enxofre 
para reduzir o volume de escórias (250- 
-275 kg esc./ton). 


2 — Temperaturas de vento mais elevadas 
O Fukuyama trabalha a 1.150ºC em compa- 
ração com os 800-1000ºC usuais noutras 
siderurgias. 


3 — Inspecção de combustível 

No Fukuyama utiliza-se 28 a 29 kg de 
óleo tonelada para reduzir o consumo de 
coque altamente dispendioso. 

(n.b. A quantidade óptima de coque com 
injecção de combustível é de 365 kgiton 
com 155 kg/ton de óleos pesados e 27º/o 
de oxigénio no vento; a melhor proporção 
sem injecção de combustível é 465 kg/ton). 


4 — Alta pressão na goela. 


O Fukuyama trabalha a 1,8-0,5 kg/cm” e 
uma pressão de vento de 2,6 kg/cm*. 


(Na U.R.5.S. usa-se mais alta pressão para 
permitir maiores quantidades de vento sem 
o perigo de criar chaminés na carga). 


5 — Controle térmico do forno por: 


a) volume de vento (resposta rápida) 
b) adição de vapor de água ao vento 
(resposta lenta) 


Fig. 4 — Encargos na produção de ferro (segundo Le- 
ckie) (7). Apresenta-se o excesso de custos de funcio- 


namento e custos fixos relativos a instalações A.F. 


e sinterização, trabalhando em condições 
de capacidade óptima, 
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c) adição de oxigénio para permitir o 
uso de mais óleo. 


Um indicador denominado Índice de Produção 
do A.F. (B.QO.I.) dado por: 


P(B + 10) 
72(D — 10) 


B.O.1I. = 


P é a produção em ton/dia, B o peso 
da carga em cwt/ton ferro e D o 
diâmetro do cadinho em pés, 


em que: 


tem sido usado para comparar a produção de 
várias unidades, quando compensadas na dife- 
rença de dimensões e tonelagem de carga. O 
Fukuyama atinge um índice 128 enquanto em 
muitos casos se considera já bom um valor de 
70 a 80. Nesta base, parece valer a pena inves- 
tigar, variando as condições operacionais, para 
aumentar a produtividade antes de considerar 
um posterior investimento de capital. Poder-se-á 
talvez aumentar de 50%/y uma actual produção 
com um dispêndio consideravelmente menor que 
o exigido por uma nova instalação. De facto 
entre 1950 e 1962 as instalações NKK Kawasaky 
aumentaram a produção de cada um dos seus 
fornos de 6,9 m de 600 para 1300 ton/dia. 


Nova aparelhagem açessória poderá ser também 
um investimento rentável, como por exemplo, 
para: 


a) análise dos gases da goela 

b) medição da temperatura das paredes 

c) medição da pressão para avaliação da 
permeabilidade 

d) análise do «stockline profile» 

e) análise contínua da carga 

f) controle por computadores para maximizar 
a produção, reduzir o consumo de coque e 
melhorar a qualidade do produto. 


2.2 — Redução directa 


Os processos de redução directa merecem 
análise cuidada pois poderá em breve ser possível 
utilizá-los econômicamente para mais baixos níveis 
de produção que os dos processos convencionais. 
Dailey (5) sugere mesmo que os custos fixos de 
produção combinada de ferro e aço podem ser 
30 a 40º/, mais baixos do que para o alto forno 
e via B.O,F., e os custos de produção inferiores 
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em 5º/o para uma capacidade de 3000000 de 
toneladas. 

Será possível usar ferro de redução directa 
de 3 modos: 


1 — Como carga de alto forno, para aumentar 
a produção, particularmente quando outras 
unidades estão em reparação ou há uma 
pressão temporária da procura. 


2 — Como carga para fornos de arco. Nas dis- 
cussões do centenário do «Iron and Steel 
Institute», * em Londres, Mr. Ascough da 
B.S.C. Colvilles (6) sugeriu que, com a 
energia a 0,0204 esc./kwh e o ferro de 
redução directa a 700-800 esc./ton, esta 
aplicação seria econômicamente atraente 
para o Reino Unido. Mr. Mekenzie da 
B.S.C. Midlands Group indicou nessa 
altura, contudo, que na sua opinião um 
preço de 2050 escudos seria mais realista. 
Obviamente a discussão âcerca das pers- 
pectivas económicas do processo continuará 
por longo tempo mas podemos já referen- 
ciar que os tempos entre os vazamentos 
podem ser considerâvelmente reduzidos por 
adopção de uma alimentação contínua a 
«pellets» (7). Isto porque se podem usar 
maiores níveis energéticos utilizando uma 
camada de «pellets» para proteger os 
refractários (fig. 5). 


3— À terceira mais importante aplicação para 
o ferro de redução directa pode ser na 
forma de pó. Haveria reduções muito 
significativas na energia necessária para a 
produção de banda se, tanto: 


a) o sistema de redução, sinterização e 
laminagem da Fundação de Pesquisas de 
Ontário (8), como: 

b) um sistema de laminagem directa do 
pó, 

se tornarem econômicamente viáveis. Isto 
depende em grande parte de uma redução 
muito considerável do custo do pó de ferro. 


* Neste encontro foi apresentado uma série de exce- 
lentes artigos de revisão que cobriam todo o campo de 
Tecnologia do ferro e aço, razão por que lhe é feita fre- 
quente referência. Estes artigos estão publicados no exem- 
plar de Junho de 1969 do «Journal of the Iron and Steel 
Institute”, 
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Uma aplicação mais imediata do pó de 
ferro de redução directa pode ser a pro- 
dução de esboços forjados, requerendo ape- 
nas um limitado grau de consolidação na 
forjagem. 


Toda a carga fundida 
a 


N 


/ 
| Periodo ce 
| refinação 
1 


so ca cargaiTfor 


Carga q 
refinação continua 


“empo à partir do cmcio do aguecimento(lhoras) 


Fig. 5 — Comparação dos tempos entre sangrias na 
fusão em forno de arco para cargas de sucata e esponja 
(segundo Sibakin, Hookings e Roeder) (7). 


3. ACIARIA 


Actualmente os fornos de soleira básicos fabri- 
cam ainda a maior parte da produção mundial 
de aço, mas esta percentagem já baixou do 
máximo de 80º/o para 55º/o, enquanto a produ- 
ção pelo processo L.D. cresceu para 25º, e a 
respeitante aos fornos de arco para 15º/,. Pro- 
vavelmente estas tendências manter-se-ão e assim 
dedicaremos a nossa atenção aos processos 
básicos usando oxigénio e fornos de arco. 


3.1 — Aciaria B. 0. F. 


Este processo desenvolveu-se muito rápida- 
mente, especialmente no Japão, onde 2/3 da pro- 
dução total de aço são presentemente obtidos por 
esta técnica. A dimensão dos fornos tem vindo 
a aumentar, de modo que temos hoje fornos de 
300 toneladas de capacidade com tempos entre 
Vazamentos inferiores a 50 minutos. Isto signi- 
fica que um conjunto de três fornos B. O.F. pode 
produzir cerca de 4,5 milhões de toneladas de 
aço por ano. À percentagem de sucata carregável 
no convertidor L.D. está geralmente entre 17 e 
30"/o. Esta quantidade é inferior às disponibilidades 
de sucata em muitos países desenvolvidos tais como 
U.S5.A. (44º/,) e U.K. (47º/,), levantando-se a 
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necessidade de usar fornos eléctricos para fundir a 
sucata disponível, pois esta é relativamente barata 
quando comparada com o ferro obtido no alto 
forno. A figura 6, obtida de um artigo de Mar- 


149 mi pg 


Fig. 6 — Custos da produção de aço em função da quan- 
tidade de metal fundido usado (segundo Marshall) (9). 


shall (9), põe em evidência o modo como os 
custos de produção do aço dependem da percen- 
tagem do metal fundido usado e das dimensões 
das unidades. 

Tal como na produção de ferro, uma das preo- 
cupações actuais é obter um melhor controle do 
processo B. O. F.. Este objectivo pode alcançar-se 
por um aumento de automatização do processo 
usando informações obtidas por monitoração do 
volume de oxigénio fornecido, temperatura, 
volume e análise dos gases de escape. Os instru- 
mentos de medida precisam ainda de ser melho- 
rados para uma medição contínua da temperatura 
e do conteúdo em oxigénio do aço. 

Estão a desenvolver-se métodos para a produ- 
ção de aços de carbono médio usando técnicas 
de descarbonização. Estas dependem da obtenção 
de curvas reprodutíveis de descarbonização 
e Meyer (10) mostrou que isso é possível se pre- 
viamente se obtiver um equilíbrio entre a oxida- 
ção e descarbonização do metal na emulsão escó- 
ria-metal-gás. No mesmo artigo Meyer mostra 
que é possível produzir aço inoxidável num con- 
vertidor L.D., embora a produção potencial 
do método seja ainda muito grande para a pro- 
cura. 


3.2 — Fornos de arco 


Estas unidades estão também a ser construídas 
com diâmetros muito maiores, até 7,2 m, e capa- 
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cidades de 180 toneladas. O rendimento destes 
fornos é cerca de 90 '/, e o tempo de operação 
inferior a 4 h (10). Por exemplo, o conjunto 
em elaboração em Templeborough, de Steel Peech 
and Tozer, tem 4 unidades de 135 toneladas cada 
e 2 de 180, produzindo 1,5 milhões de toneladas 
anuais de aço com um consumo de energia 
de 520 Kwh por tonelada. 

Para obter produções elevadas deve fazer-se 
o possível para minimizar as demoras. Isto 
envolve: 


1 —um controle muito cuidadoso do enfarda- 
mento de sucatas e carga; 


2 — programação de potência a fornecer condi- 
cionada pelas restrições máximas da pro- 
cura ; 


3 — utilização do oxigénio do modo mais efec- 
tivo ; 

4 — melhoramento dos materiais para eléctro- 
dos e das técnicas de manuseamento ; 


5 — cuidado na especificação e manutenção dos 
refractários, usando um programa estudado. 


As tendências futuras orientam-se no sentido 
da utilização de fornos de mais alta potência, 
usando mais elevadas razões de transformação, 
seja, 75000 — 100000 KVA, que devem dar 
origem a maiores velocidades de fusão, até 
70 ton/h. Se os preços para a obtenção de «pel- 
lets» se tornarem razoáveis, serão provavelmente 
utilizados como carga continua para aumentar 
a velocidade de produção. Foi sugerido (11) que 
uma carga de 60 "/, de «pellets» e 40 “/ de sucata 
seria uma combinação razoável, na medida 
em que é necessário usar sucata proveniente 
de fábricas e sucata proveniente do mercado 
a preço razoável. O futuro dos fornos de arco 
parece bom, particularmente no que diz respeito 
à possibilidade de obter energia mais barata 
de origem nuclear. Se esta previsão não for 
totalmente realizável, alguns produtores consi- 
derarão possivelmente a hipótese de usar quei- 
madores a gás ou Óleo no período de fusão. 
Os custos de produção para um processo SLRN 
e forno de arco foram estimados em 870 esc./ton 
para uma produção anual de 1,5 milhões de tone- 
ladas. 
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3.3 — Novos processos 


Uma série de processos de aciaria contínuos 
tais como: 1) «Spray Steelmaking process» 
da BISRA, 2) o processo IRSID ou 3) o pro- 
cesso WORNER que está a ser experimentado 
em Rio Tinto, estão presentemente a ser apre- 
ciados. 

Todos estes processos se baseiam no aumento 
da área de contacto entre o banho metálico 
e a escória em presença do oxigénio, o que se pro- 
cura conseguir de diferentes modos. No pro- 
cesso BISRA desintegra-se uma corrente de metal 
líquido por jactos de oxigénio e cal. No pro- 
cesso IRSID injecta-se oxigénio na superfície 
do banho metálico contido num recipiente 
e obriga-se a emulsão escória-banho a escorrer 
sobre um dique para um recipiente de decanta- 
ção. O processo WORNER utiliza um tanque 
comprido onde o banho e a escória circulam 
em contra corrente. Fazem-se incidir jactos 
de oxigénio sobre a superfície de modo a formar 
uma emulsão na parte do forno que constitui 
o reactor própriamente dito. 

Muito há ainda para fazer no desenvolvimento 
de métodos adequados de controle para estes 
processos. Num futuro ainda distante parece 
viável que o processo contínuo será usado para 
o grosso da produção de aço, mas não parece ter 
suficiente flexibilidade para o pequeno produtor. 
Num futuro mais longínquo deve dar-se atenção 
à possibilidade de usar energia nuclear na 
produção de aço. O Prof. Antonio Antonioli (6) 
descreveu um projecto de estudo em que a E.E.C. 
está envolvida, para uma instalação usando 1/3 
do calor fornecido por um reactor nuclear refri- 
gerado a gás para uma redução directa, possivel- 
mente num reactor de leito fluidizado, usando 
gás natural, sendo os restantes 2/3 para produzir 
energia eléctrica destinada a fornos de arco para 
produção de aço a partir do ferro de redução 
directa. 


4. PROCESSAMENTO DO AÇO 


As duas modificações mais importantes neste 
campo tem sido o desenvolvimento de unidades 
de grandes dimensões para desgasificação em 
vazio e vazamento contínuo. 
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4.1 — Desgasificação em vazio 


EA 


A desgasificação em vazio é presentemente 
mais significativa na produção de aços ligados 
para remoção de hidrogénio em grandes lingotes 
e controle de inclusões (fig. 7 e 8) (12). 


Fig. 7 — Desgasificação em vazio no vazamento 
para o molde. 


De acordo com Holden (13) os custos são da 
ordem de 28 a 70 esc/ton. Alguns produtores 
estão a começar a encarar as técnicas DHHU 


Pressao mm bg 


Hidrogênio mi Hoogm 


Fig. 8 — Conteúdo em hidrogénio obtido durante a des- 
gasificação em vazio (segundo Speith, von Ende 
e Pfeiffer) (12). 
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e RH (fig. 9 e 10) como fornecendo um novo 
estágio ma desgasificação e ajustamento da 
composição na produção do aço (fig. 11 e 12) 
(14 e 15). Este novo processo parece ter impor- 


Desgasificação DH 


Tremonha dos 
elementos de 


liga vazio 


Elemento de aquecimento 


Camara de vazio 


Fig. go — Representação esquemática de uma unidade 
de desgasificação D.H. 


tância para o grosso da produção de aço, no 
desenvolvimento de substitutos de aço eferves- 
cente para vazamento continuo. A Companhia 


Desgasific ação RH 


Adição de elementos de liga 
Vazio 


Elemento 
de ” 
aquecimento 


Camara de 
vazio 

+ *s 
Saida de gas 


Fig. 10 — Representação esquemática de uma unidade 
de desgasificação R.ll. 
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Fig. 11 — Desvio das condições de equilíbrio na reacção 
carbono-oxigénio verificadas em vários processos de 
aciaria (segundo Pearson, Venkatari e O Hanlon) (14). 


de Aço de Wales está já a instalar uma unidade 
DH que terá uma capacidade de 1000000 
ton/ano. 
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Fig. 12 — Desvio das condições de equilíbrio na reacção 
carbono-cxigéênio nas operações de desgasificação em 
vazio (segundo Sieckbert) (15). 


4.2 — Vazamento contínuo 


Há economias evidentes no vazamento con- 
tínuo por causa do aumento de produção e 
eliminação dos fornos «pits» e trens de lami- 
nagem primários. Presentemente existem cerca 
de 300 unidades na indústria do aço, algumas 
das quais permitem manusear 200 toneladas de 
aço de cada vez. Algumas unidades trabalham 
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Os Mobiltac 


aguentam-se . 
tanto em servico, 


que dificilmente 
poderemos . 
vendê-los muito. 


Os Mobiliac constituem uma linha de 
produtos criada especificamente para lubri- 
ficação de engrenagens abertas. Por oulras 
palavras, o desperdício de Óleo é mínimo 
Ou mesmo inexistente. 

Isto significa menor frequência na apli- 
cação, maior protecção das engrenagens 
e - infelizmente para nós - menos consumo 
de lubrificante. 

Mas a grande adesividade não é a única 
vantagem dos Mobiltac. Outra, é o seu 
bom funcionamento a baixa temperatura 
mesmo abaixo de zero continua a prole- 
ger eficazmente as engrenagens, evitando 
avarias. 

a E Peça à MOBIL mais informações sobre 
on O SI NUR O Cn estes produtos e poderá economizar anual- 
dc a mente milhares de escudos. 


Mobil 


vendemos mais vendendo menos 
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A Mobil 


mostra-lhe 


inumeras maneiras 
de gastar 
menos óleo. 


a 


Por ano poupamos à 
industria varios milhares 
de escudos. 

Uma grande parte 

deste dinheiro 

sai das nossas algibeiras. 
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Os Mobiltac 
ntam-se 
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que dificilmente 

poderemos 


vendé-los muito. = 


vendemos mais vendendo menos 


vendemos mais vendendo menos 
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Este simples ensaio 
mostra como reduzir 

os custos 

de manutenção 

dos sistemas hidráulicos. 
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o seu dinheiro! 


vendemos mais vendendo menos 


Se é raro apresentarmos 
o preço mais baixo, 
porque razão vendemos 
tanto óleo para 

grandes turbinas? 


vendemos mais vendendo menos 


«slabs» com secção máxima de 2600 cm”, mas 
três companhias podem obter por este processo 
secções até 6500 cm?. As velocidades de vaza- 
mento variam entre 0,5 e 6 m por minuto para 
os «billets» mais pequenos. Em termos de tone- 
lagem isto quer dizer que são possíveis cerca de 
20 ton/h para cada «strand» com billets de 
250 mm e 150 ton/h para os «slab» de maiores 
dimensões. As produções anuais variam de 
10000 ton/ano para uma unidade de pequena 
secção apenas com um «strand», até 800 000 
ton/ano para unidades com 8 «strand» e grande 
secção. 

Técnicas envolvendo a utilização de moldes cur- 
vos ou flexão e rectificação do «strand» reduziram 
a altura deste equipamento (fig. 13) e os melho- 


Molde curvo” 


Rolos quias 


Fig. 13 — Representação esquemática do vazamento 
contínuo usando a técnica do encurvamento do «strand». 


ramentos futuros parecem estar nas técnicas de 
redução do «strand», moldes de largura variável 
e sequência de vazamento. Nas técnicas de 
redução do «strand», os rolos guias comprimem 
o material vazado quando neste resta um núcleo 
líquido de 3 a 5% da sua secção. Esta acção 
obriga o metal fundido a voltar ao molde e 
parece reduzir substancialmente a porosidade 
e segregação centrais. Uma posterior redução de 
secção aumentará a velocidade de produção por 
«strand» para pequenos «billets» e permite a 
obtenção de variados tamanhos finais para cada 
tamanho de molde. 
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A redução do «strand» está já em aplicação 
na Badische Stahlwerk onde a produção de dois 
fornos de arco de 50 toneladas cada é vazada 
numa unidade de 6 «strand». As dimensões 
iniciais são 90>< 130 mm posteriormente redu- 
zidas para 90><90mm antes da laminagem 
em varão para betão armado. 

A firma Algona vazou continuamente 620 
toneladas através de 4 «strand» em 6h 30 min., 
e esta tendência manifesta-se ainda numa 
reivindicação russa de que se conseguiu vazar 
um «slab» com 1 km de comprimento. 

Os custos iniciais citados por Holden (13) 
foram avaliados em 350-650 esc/ton/ano e os 
custos de operação são cerca de 70 a 140 esc. 
mais baixos que os referentes ao processo da 
produção de lingotes com utilização de trem de 
laminagem primário. 

A grande vantagem do vazamento contínuo é 
o corte da linha de produção. Reduções por 
laminagem de 6:1 até 10:1 consolidam a 
estrutura central e dão origem a propriedades 
similares às obtidas em lingotes tratados conven- 
cionalmente com reduções de 10:1 a 250:1. 
Esta vantagem verifica-se também para «blooms» 
perfilados e a firma Algona, cooperando com a 
BISRA, já teve em funcionamento durante um 
ano inteiro uma unidade de 2 «strand» para 
vazamento de «blooms» perfilados. Os alonga- 
mentos na laminagem são de 9:12 12:1 face 
aos 55 a 57 :1 verificados nos processos conven- 
cionais. 

Os maiores problemas que limitam as veloci- 
dades de produção residem principalmente na 
resistência dos materiais refractários utilizados 
nos «buzettes» e panelas de vazamento. 


5. TRABALHO MECÂNICO 


A tendência para o corte das linhas de produ- 
ção exemplificada pela utilização de vazamento 
contínuo está também patente na atenção posta 
na laminagem de pós. 

Embora a aplicação desta técnica a aços de 
baixo carbono dependa de um abastecimento 
barato de pó de ferro, há uma perspectiva mais 
imediata na aplicação do processo à produção de 
banda de aço inoxidável. A laminagem da banda 
em aço inox é um processo trabalhoso visto que 
as reduções que se podem efectuar após recozi- 
mento são pequenas. 
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Portanto talvez se possa pagar pó de aço inox 
produzido por um processo de atomização e ter 
ainda possibilidades de produzir banda de aço 
inox economicamente. 

Aparte estas mudanças radicais nos processos 
de trabalho mecânico, a tendência actual reside 
em aumentar a utilização racional das instalações 
e reduzir as perdas e rejeições através de um 
controle automático de dimensões. Isto evidente- 
mente significa melhoramento na produção porque 
para volumes actuais de produção muitos trens 
de laminagem operarão a 40-50º/y do seu tempo 
máximo de serviço. 

O desenvolvimento técnico no que diz respeito 
à laminagem de varão e barra orienta-se para 
processos usando «strand» múltiplos, trens mais 
rígidos e maior controle automático. Por exemplce, 
Green (16) sugeriu um trem de laminagem para 
varão com 4 «strand» produzindo varão com 5 
a 12,5 mm de diâmetro. As técnicas de «strand» 
múltiplos são usadas para as etapas grosseiras e 
intermédias e nos estádios finais usam-se 4 
«strand» simples com velocidades de acabamento 
até 3 m minuto e uma produção até 500 000 
ton/ano. Os «billets» iniciais têm até 100 mm? 
de secção e um comprimento de 12-18 m. 

Juntamente com estes melhoramentos temos 
modificações nos fornos incluindo a introdução 


Unidades 
5 altos fornos (475 m) e unidades 
associadas de sinterização e manu- 


seamento do minério . 


0 convertidores L. D. (200 ton.) . 6,4 
2 fornos de arco de 90 ton | 0.36 
2 fornos de arco de 0 ton . | 1,5 
Vazamento contínuo (Holden) (13) . | 
Trem universal de «slabs» +. +... +53 
Trem simples de chapa 210 cm .... 0,3 
Trem simples de chapa 330 em . .. | 0,3 
Trem para banda a quente (205 em) . 3,25 
Trem para banda a frio (205 em) .. (1,8 
Trem de «blooms» (107 em) .... 1,55 
» é ADO sas == sau LO 
» a fEMMsrcissm,a.G 1 0,25 
» à VAO «essnvymsna ns 0,25 
» » Secções médias ..... | 0,4 


Produção | 
(milhões de ton/ano) | (milhões de escudos) | 


do «walking beam» e «notched hearth» que são 
muito mais flexíveis que os «pusher furnaces» 
quando há atrasos na laminagem. 

Nas últimas operações da produção estão a 
ser usados «open coil laying» e arrefecimento 
controlado ao ar, para melhorar a estrutura dos 
produtos para aplicações em trefilagem. A 
tendência para as grandes unidades é mais marcada 
nos trens de laminagem de banda, a quente, 
produzindo 3 a 4 milhões de toneladas por ano. 
Isto trouxe a necessidade de aumentar a velocidade 
de decapagem, usando, por exemplo, ácido 
clorídrico em vez de sulfúrico, e, para a chapa 
estanhada, a necessidade de prover unidades 
continuas de recozimento e revestimento. 

Estas unidades podem trabalhar agora até 
600 m'minuto. Os custos iniciais destas unidades 
são elevadíssimos e não podem portanto ser 
subutilizadas. O quadro da fig. 14 mostra os 
custos aproximados das unidades de produção de 
aço e ferro segundo Lekie e Morris (2) e põe em 
evidência os altos custos dos trens de laminagem 
a quente para a produção de banda. 


6. PRODUTOS 


Os progressos na metalurgia dos materiais fer- 
rosos continuam com grande expansão, mas men- 


Custos | Custos por ton/ano 


(escudos/ton) 


| | 


5.250 1.050 
3.220 500 


200 TIO 
700 | 470 
350-560 

USO 218 
780 2.600 
1.099 2.198 
2.300 TIO 
030 790 
350 | 226 
560 | 500 
420 1.680 
350 1.400 

| 510 | 1.270 


Fig. 14 — Custo de várias unidades siderúrgicas (segundo Leckie e Morris) (2) 
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cionaremos apenas três campos onde o seu efeito 
na tecnologia da produção é susceptível de ser 
mais significativo. Estes são: 


1 — Aplicação de ciclos rápidos de recozimento 
para obter grão fino em banda, ex. 6 8 
ao recozer, após 85 a 90º/o de redução a 
frio (aquecendo a 700ºC em - 30s e arrefe- 
cendo de novo até 200 em 1005.). Esta 
técnica leva-nos à obtenção de uma tensão 
de cedência de 41 kg/mm” na folha esta- 
nhada (0,07º/, €.) 


2 — Utilização de laminagem controlada para 
produzir aços estruturais de alta resistên- 
cia. Nesta técnica controla-se de novo o 
tamanho do grão. Neste caso a resistência 
depende de dois factores: primeiramente 
adicionando 0,03º/, de nióbio obtém-se 
carbonetos de nióbio estáveis que ajudam 
a controlar o tamanho do grão; depois um 
controle cuidadoso dos esquemas de aque- 
cimento e laminagem dá azo a que certa 
fracção de carboneto de nióbio se dissolva 
a alta temperatura e se precipite na matriz 
ferrítica de grão fino (fig. 15) (17). Estes 
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Fig. 15 — Tensão de cedência dos aços contendo Nió- 
bio em função do tamanho de grão e tratamento tér- 
mico (segundo Irvine) (17). 


aços estruturais de alta resistência com 


tensões de cedência até 50 kg/mmº pode- 
rão ser obtidos a preços para os quais a 
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economia em peso numa estrutura facil- 
mente compensará o excesso de custo 


requerido se os trens forem concebidos 
para suportar os novos esquemas de lami- 
nagem. Actualmente a interrupção dos 
ciclos normais dá azo a encargos elevados. 


3— A utilização de técnicas de extrusão e forja 
a frio para enformar aço com 0,1º/, €. 
Esta técnica economiza muito material 
quanto comparada com a maquinação e é 
a melhor para enformar um grande número 
de componentes industriais. 


7. CONCLUSÕES 


Em conclusão: — Está a processar-se uma 
rápida evolução na indústria do aço e a procura 
para o material parece continuar a aumentar. O 
controle é de importância vital na administração 
e na tecnologia. Um país como Portugal, que 
expandirá provavelmente a sua indústria siderúr- 
gica na próxima década, terá muitas e fascinan- 
tes decisões a tomar. 


As mais importantes serão: 


1— Até que ponto o ferro de redução directa 
e os fornos de arco fornecem uma via eco- 
nómica de produção de aço. 


2— Até que ponto o vazamento contínuo 
substituirá os processos convencionais uti- 
lizando o lingote e a laminagem primária 
para obter «blooms» e «slabs», 


3 — Até que ponto o controle dos processos de 
laminagem será usado para obter proprie- 
dades préviamente conseguidas através de 
tratamentos térmicos ou elementos de 
liga. 
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TEORIA DA QUANTIFICAÇÃO DAS ACTIVIDADES 
PROFISSIONAIS E SUA APLICAÇÃO 
À ENGENHARIA ELECTROTÉCNICA 


RESUMO 


Atraves da noção de plenitude, que traduz a relação 
entre uma actividade e os problemas gerais da respectiva 
profissão, estabelece-se uma teoria de quantificação, que 
aplicada à engenharia electrotécnica ordena as actividades 
disponíveis na sociedade. 


1. INTRODUÇÃO 


Um engenheiro recém-formado, após a forma- 
ção escolar, encontra-se numa situação de escolha 
entre várias actividades a que se pode dedicar. 
Pela carência de elementos no julgamento sobre 
a direcção a tomar parte da situação embaraçosa 
para uma resolução baseado frequentemente em 
factores secundários, decisão por vezes com 
reflexos desagradáveis, mais tarde, numa altura 
em que já se torna difícil a mudança de rumo. 

Procura-se neste estudo comparar quantitati- 
vamente as diversas actividades possíveis do 
engenheiro electrotécnico na nossa estrutura 
actual segundo um critério definido a partir da 
maior ou menor participação nos problemas gerais 
da profissão. A valorização deve pois ser encarada 
sob este aspecto. Tecem-se depois considerações 
sociais a partir da teoria estabelecida. 


2. DOMÍNIO DAS ACTIVIDADES PROFIS- 
SIONAIS 


A noção de actividade profissional corresponde à 
ocupação do homem, em relação à sua profissão, 
na sociedade em que está integrado. 

O domínio de actividades profissionais inerente 
a uma dada profissão é mais ou menos extenso 
conforme engloba mais ou menos actividades 
possíveis nitidamente bem caracterizadas e, em 
certa medida, bem diferenciadas. Portanto o 
domínio das actividades profissionais de uma profis- 
são é definido pelo conjunto ordenado, segundo 
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SYNOPSIS 


The relation between an activity and the essential 
problems of a profession makes possible the quantification 
of the profession. By means of thns theory the electrical 
engineering activities are ordened, 


um determinado critério de hierarquia, das activi- 
dades profissionais relativas a essa profissão. 

Introduz-se uma abstracção na teoria, justifi- 
cada pela generalidade e simplicidade do trata- 
mento matemático, considerando o domínio pro- 
fissional, caracterizado entre certos limites, como 
um conjunto contínuo de actividades. Esta 
abstracção vai permitir tratar este conjunto como 
continuo e por consequência representável gráfi- 
camente sobre uma recta. 

Na prática todavia apenas existe um número 
discreto de pontos nitidamente diferenciados rela- 
tivos às actividades disponíveis na sociedade. 
Admite-se então o 1.0 postulado da teoria: as 
actividades profissionais diferenciadas, estando 
ordenadas entre os limites do domínio profissio- 
nal segundo um determinado critério de hierar- 
quia, estão igualmente afastadas entre si. Deste 
modo se quantificam as actividades entre si. 

Assim se a uma profissão corresponde um 
domínio profissional entre a e b no domínio 
genérico das actividades profissionais é possível 
distribuir as n actividades diferenciadas disponí- 
veis na sociedade inerentes a essa profissão no 
intervalo (a, b), determinando-se pontos igual- 
mente repartidos nesse domínio rectilíneo e limi- 
tado (fig. 1). 


A À b 
peca e ha E» ba AS 
0 1 2 3 n- n 
À —s 


Fig. 1 — Representação gráfica linear do domínio de 
actividades profissionais de uma profissão, 
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Exteriormente ao intervalo (a, b) não há acti- 
vidades ligadas à profissão. Então pode escolher- 
-se, em qualquer caso, o limite a para origem 
de ordenação das actividades e atribuir-lhe o 
valor nulo. 

Em toda a profissão admite-se que existe uma 
actividade perfeita, correspondente a uma noção 
arquétipa, com a qual o homem completamente 
formado dentro da sua profissão teria oportuni- 
dade de trabalhar com todos os problemas, de 
qualquer natureza, científica, ética, etc., que a 
profissão suscitasse. Esta actividade perfeita cons- 
tituirá um limite que se toma para origem da 
distribuição de actividades no domínio da pro- 
fissão e que será atingida pela noção limite de 
perfeição (dentro do conceito exposto), portanto 
transcendente. 

Depois existem actividades reais, sucessiva- 
mente mais imperfeitas por se referirem a aspectos 
cada vez mais restritos da profissão no que diz 
respeito aos problemas fundamentais dessa pro- 
fissão. 

Há-de existir uma actividade total, caracterizada 
pela presença de todos os problemas fundamen- 
tais ligados à sociedade sem os transcender. 

E, por fim, surgirá a actividade de subsistência a 
partir da qual o homem já não exerce actividade 
dentro da sua profissão, constituindo aquela acti- 
vidade a última que lhe permite ainda conside- 
rar-se na sua profissão, mas já extremamente 
rarefeita. 


3. PLENITUDE PROFISSIONAL 


Para ordenação das actividades dentro do seu 
domínio é necessário definir um critério de valo- 
rização que permita quantificar as diversas acti- 
vidades de uma profissão. 

Ora uma profissão é caracterizada por um certo 
número de problemas fundamentais que dão forma 
à profissão e que estão na base do conjunto das 
actividades. 

Cada uma das actividades pode estar directa- 
mente ligada a um ou mais dos problemas fun- 
damentais, pode encontrar-se relacionada directa- 
mente com um, menos acentuadamente com outro 
e desligada dos restantes, ou pode haver, enfim, 
qualquer outra combinação possível com o número 
de problemas gerais. 

Para se quantificarem as actividades de uma 
profissão admite-se o 2.º postulado da teoria que 
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atribui a cada actividade um valor, numa escala 
arbitrária, conforme se ocupa de cada um dos 
problemas fundamentais mais ou menos directa- 
mente. 

Deste modo estabelece-se uma matriz de coefi- 
cientes de actividade entre cada actividade e cada 
problema fundamental (Quadro 1): o coeficiente 
C.. traduz a ocupação directa da actividade j com 
o problema k. Este número é absoluto mas só tem 
significado relativamente aos restantes, dada 
a arbitrariedade da escala de valores. 


Tipos de nn “Problemas fundamentais a 
activ. | B, B ... Bo Totais 
A | Cy Gm um Kyg 1 
Ao Ca Cy Cop C 
An Ca Cn? Cnp Cn 

Totais | CG CG ... CG | € 


QUADRO I — Coeficientes de actividade. 


A actividade genérica A, terá portanto um peso 
total igual ao somatório dos coeficientes de acti- 
vidade relativos a essa actividade e ligados 
a todos os problemas gerais 


A (1) 


k==] 


Efectuando a soma total dos pesos de todas 
as actividades obtém-se o peso da profissão 
na escala arbitrada 


C=2€c (2) 


que se pode fazer igual a 100 */o. Por este pro- 
cesso determinam-se os valores relativos dos 
pesos das actividades em relação ao peso da pro- 
fissão, expressos em percentagem 


C, 
T, = -100 (3) 


que se denomina plenitude da actividade profissional 
pois traduz a relação directa de um profissional 
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com a totalidade dos problemas fundamentais 
da sua profissão. 

Por este meio retirou-se toda a influência dos 
valores absolutos arbitrados para os coeficientes 
de actividade permanecendo tão sômente o signi- 
ficado relativo de uns valores para os outros. Por 
isso se exige, para a conveniente aferição de umas 
pesagens em relação às outras, o profundo conhe- 
cimento dos problemas fundamentais da profis- 
são, além das diferentes actividades disponíveis 
na sociedade. Este é pois um aspecto melindroso, 
susceptível de controvérsia e exigindo necessária- 
mente permanente crítica na evolução temporal. 

Importa notar que o conceito de plenitude nada 
tem que ver com o de felicidade, pois a quanti- 
ficação introduzida refere-se apenas à maior 
ou menor presença dos problemas fundamentais 
da profissão em cada actividade dessa profissão, 
mas não se encontra muito desligado do senti- 
mento de «realização profissional». 


4. EQUAÇÃO DA PROFISSÃO 


Estabelecidas as bases anteriores podem cons- 
truir-se vários desenvolvimentos. Não nos inte- 
ressam, por ora, dado o fim em vista ser a deter- 
minação de uma lei de correlação entre as activi- 
dades e as plenitudes. 

Salienta-se porém que a quantificação postu- 
lada permite estudos de vária natureza utilizando 
a análise matemática. 

Cita-se, por exemplo e porque se trata de uma 
generalização unificadora, a dedução da equação 
de uma profissão. 

Para tal basta referir o quadro matricial ante- 
rior mas com os coeficientes de actividade redu- 
zidos aos valores porcentuais, ou seja às pleni- 
tudes (Quadro II). 


Activi- Problemas fundamentais 
ci B, B cc. B Totais 
A| Hj Fo Tip É 
An | Ta Too Tp da 
Totais T, E, E | 


QUADRO II — Coeficientes de plenitude. 
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Definindo o vector actividades 


A= |. (4) 


e o vector problemas fundamentais 
B=[T,T,;...T,| (5) 


com a matriz dos coeficientes de amplitude como 
matriz de transferência 


Ta Ego vo Tg 
E Tg «Tê 

T = o . o (6) 
Tas é : d Ega 


deduz-se a equação matricial da profissão 
B==T.A (7) 


corréspondente ao sistema de equações da profissão 
que relacionam cada problema fundamental com 
as diversas actividades por meio de uma combi- 
nação linear cujos coeficientes são as plenitudes 
da matriz de transferência 


B, ca T A, + Ty A, ora Tn A, (8) 


onde k=1,2,. . .,p. 

Como os problemas gerais se quantificaram 
pelos coeficientes T, e as actividades pelos coefi- 
cientes T, será igualmente 
com k=1,2,. . .,p. 

A equação (7) condensa de um modo simpli- 
cissimo a questão complexa da quantificação da 
profissão. 


5. FUNÇÃO PLENITUDE - ACTIVIDADE 
Tomando a plenitude T como uma variável 


contínua e análogamente para a actividade A, 
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a função plenitude-actividade define-se num 
sistema de eixos rectangulares em que o eixo das 
abscissas constitui o domínio das actividades 
profissionais e o eixo das ordenadas as pleni- 
tudes profissionais respectivas. 

As Pplenitudes, sendo valores quantificados, 
estão ordenadas na parte positiva do eixo. 

Falta ordenar as actividades dentro do respec- 
tivo domínio na parte positiva do eixo das 
abscissas. Para isso estabelece-se o critério de 
hierarquia das actividades: o ordenamento das 
actividades de uma profissão corresponde ao 
decréscimo da plenitude profissional. 

Em consequência deste critério a função pleni- 
tude-actividade é uma lei hiperbólica (fig. 2) 


(10) 


onde T é a variável plenitude, A a variável activi- 
dade e T, uma constante que depende da natureza 
da profissão. 


A —e 


Fig. 2 — Representação gráfica da função plenitude- 
-actividade T(A). 


Considere-se o domínio das actividades com- 
preendido entre O e A,, o domínio das pleni- 
tudes entre O e T,, designando-se A, e T, valores 
extremos da profissão. 

Define-se a função fronteira linear 


A 
Teo-L (A=8) (11) 
A 


o 
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que intersecta os eixos coordenados em A, eT,. 
Esta recta limita a zona do plano plenitude-acti- 
vidade em que é possível haver correspondência 
entre as actividades e as plenitudes, para uma 
profissão cujos limites de actividades e de pleni- 
tudes sejam respectivamente A, e T,. 

As intersecções das curvas expressas por (10) 
e (11) determinam os pontos concretos limites do 
espaço do plano em que correspondem activida- 
des concretas, isto é, com significado científico 
e social. 

Ao ponto concreto limite relativo à maior ple- 
nitude chama-se ponto total por corresponder à 
actividade total, quer dizer, com uma plenitude 
de 100º/,. 

O ponto concreto limite correspondente à 
menor plenitude poderá designar-se ponto de subsis- 
tência por corresponder à actividade de subsistên- 
cia praticável na sociedade com o mínimo de 
plenitude sem deixar de exercer a profissão. O 
valor da plenitude de subsistência é função dos 
problemas fundamentais da profissão e do número 
de actividades diferenciadas. 

O cálculo analítico destes pontos conduz à 
equação característica da profissão 


ô T 
A-AÇÃ+A mo =0 (12) 


o 


sendo A real na condição A, > 4T,/T,. No caso 
em que A, =4T,/T, só há uma solução 
(A=A,/2), isto é, a recta da função fronteira 
é tangente à hipérbole da função plenitude- 
-actividade, o que corresponde a haver somente 
ama actividade praticável para essa profissão. 
Portanto a essa actividade única correspondem 
todos os problemas fundamentais e para além 
dela muda-se de profissão. 

As actividades dentro do intervalo [0,4,] 
aberto à direita são transcendentes, definindo o 
campo transcendente de actividades, pois correspondem 
a valores de plenitude superiores a 100/,. 

As actividades contidas no intervalo [A., A,) 
aberto à esquerda do domínio de uma profissão, 
relativos a valores de plenitude inferiores ao mi- 
nimo praticável (inferiores à plenitude de subsis- 
tência), são actividades de mudança de profissão. 

Deste modo o domínio de actividades de uma 
profissão aparece constituído por três campos: 
transcendente, concreto e degenerescência (fig. 3). 
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Campo Campo Campo de 

transcendente concreto degenerescência 

a E 
A —— 


Fig. 3 — Campos do domínio de actividades, 


Os valores extremos A, e T, são os que cor- 
respondem à função fronteira. Mas nestes cam- 
pos a função T(A) excede a fronteira. Resulta 
daqui o 3.º postulado da teoria: a variação da fun- 
ção plenitude-actividade nos campos transcen- 
dente e de degenerescência é linear e de acordo 
com a função fronteira. 

Estabelecidos estes conceitos fundamentais 
considere-se uma aplicação para exemplifica- 
ção. 


6. PROBLEMAS FUNDAMENTAIS DA ENGE- 
NHARIA ELECTROTÉCNICA 


À engenharia electrotécnica apresentam-se 
os seguintes problemas essenciais: I — Concepção 
de equipamentos e de instalações; II — Fabrico 
de equipamentos; III — Instalação de equipamen- 
tos; IV Exploração e manutenção de equipamentos 
e de instalações. 

Ao engenheiro electrotécnico cumpre criar 
e resolver as implicações destes problemas, cons- 
tituindo assim a sua actividade completa dentro 
do domínio normal da profissão. 


7. DOMÍNIO DE ACTIVIDADES DO ENGE- 
NHEIRO ELECTROTÉCNICO 


São os problemas gerais definidos anterior- 
mente que dão forma à profissão do engenheiro 
electrotécnico. As actividades socialmente estabe- 
lecidas que se preocupam com esses problemas 
constituem a população discreta do domínio 
de actividades da profissão em estudo. 

O Quadro III indica os tipos de actividade 
disponíveis ao engenheiro electrotécnico, dentro 
da nossa sociedade, tendo-se estimado os coefi- 
cientes de actividade em relação a cada problema 
fundamental, que conduzem ao cálculo das res- 
pectivas plenitudes. 
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8. FUNÇÃO PLENITUDE-ACTIVIDADE DO 
ENGENHEIRO ELECTROTÉCNICO 


De acordo com esta análise conclui-se: 


a) O ensino corresponde à actividade total, pois 
a plenitude correspondente é de 100º, visto 
ocupar-se indistintamente dos quatro problemas 
gerais ; 

b) A investigação científica ocupa o lugar 
seguinte na ordem de plenitudes decrescentes, 
apesar de se ocupar essencialmente apenas de um 
dos problemas fundamentais ; 

c) A montagem de instalações aparece a seguir 
como sendo a actividade prática mais completa ; 

d) A montagem de fabrico e a exploração sur- 
gem com quase igual valorização embora se 
ocupem de problemas distintos; 

e) Os projectos de instalações e de fabrico têm 
Obviamente um valor semelhante em ocupações 
de natureza complementar ; 

f) A administração desce levemente, sendo 
dificil generalizar a correspondência com cada um 
dos problemas dada a diversidade de tipos de 
empresa pelo que se atribuiu um peso mais ele- 
vado no problema de concepção. 

g) O comércio situa-se no lugar mais reduzido. 


Estes resultados podem inferir-se imediata- 
mente após a quantificação estabelecida pela 
teoria de quantificação das profissões anterior- 
mente concebida. 

Defina-se o 2.º postulado dando à predominân- 
cia de um problema fundamental numa activi- 
dade o valor arbitrário 25, que se refere à 
percentagem visto que existem quatro problemas 
gerais. Analise-se agora cada correspondência 
actividade-problema e atribua-se um peso relati- 
vamente aos anteriores, construindo-se assim a 
matriz de plenitudes. 

A representação gráfica da função plenitude- 
-actividades baseando-se no 1.º postulado situa 
as actividades igualmente afastadas no domínio 
das actividades, e segundo o critério de hierar- 
quias coloca as actividades por ordem decrescente 
de plenitudes. 

Definiram-se nove actividades diferenciadas 
possíveis dentro da estrutura social. Introduzindo 
a meio do campo transcendente a actividade filo- 
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ACTIVIDADE 


TIPO DE SERVIÇO ONDE 


| PROBLEMAS GERAIS 


TIPO o NATUREZA SER TERNO A 6 ANA rfp | IV 
| F E , | 
Ensino 4 erelinetçd Escolas técnicas médias e superiores 25 | 25 | 25 | 25 | 100 
Investigação | Estudos e ensaios Centros de Investigação. Laboratórios || 50 125/25 0 65 
| . Secção de projecto em fábricas de ma- E E 
o Projecto de fabrico terial eléctrico 25 25] O 0 27,9 
Indústria E E | | O 
Gabrico) | Mocingoo de Linhas de fabrico em fábricas de mate- : ' 
| : : 5 2 [28 | Tô 40 
| fabrico | rial eléctrico 
o o | pet FER ane mr o PL Yi. 
; ' . Serviço técnico-comercial de empre- 
e | ppsiememacirins | sas de venda de material eléctrico. IO | 2,5 3 2,5 20 
ivenda é ici Serviço de aquisição de materiais 
o Gabinete de projectos para concurso | 
Instalação Projecto Gabinete de estudos para execução 23/25/25) 0 30 
(compra e venda)|——————— — a | td, HE, ES 
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sofia da ciência electrotécnica, a variável teórica- 
mente contínua Á (actividade) acumula-se em 
dez pontos que caracterizam as dez actividades 
diferenciadas estabelecidas. 

À actividade de subsistência A, == 10 corres- 
ponde a plenitude T,==20, pelo quadro III. 
Então a função hiperbólica, que terá de passar 
por este ponto, determina a constante T, = 200. 
Deste modo a equação da função plenitude-acti- 
vidade da engenharia electrotécnica é 


La 200 (13) 


e portanto no ponto de actividade total, para 
T,== 100, será A,==2, quer dizer a segunda 
actividade definida representa a actividade total, 
como era de esperar pois se introduziu uma 
actividade no campo transcendente. Efectuando 
a representação gráfica de função (13) cons- 
tata-se um perfeito acordo com os valores arbi- 
trados para as plenitudes (fig. 4). 
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Fig. 4 — Função plenitude-actividade da 
electrotécnica. 
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9. CONCLUSÕES 


A comparação da curva teórica, deduzida 
da teoria precedente, e os pontos relativos aos 
coeficientes estimados pela análise da natureza 
de cada actividade diferenciada e cada problema 
fundamental leva a concluir sobre a importância 
do método. 

São discutíveis as pesagens atribuídas às acti- 
vidades, pois o conhecimento dessas relações 
é em certo grau subjectivo. Contudo conjugando 
o raciocínio com um método estatístico, através 
da aplicação a uma população significativa, 
poder-se-ão determinar resultados bastante acei- 
táveis para cada profissão. 

É igualmente evidente que as conclusões deter- 
minadas numa época histórica tenderão a modi- 
ficar-se, por um lado devido à evolução da socie- 
dade, e portanto da natureza das actividades 
diferenciadas, e por outro lado devido à evolu- 
ção da profissão que com base no progresso cien- 
tífico altera os problemas fundamentais. 

Dentro dos princípios enunciados convém notar 
a ordenação das actividades que se apresentam ao 
engenheiro electrotécnico dentro do campo con- 
creto: 


a) Ensino 

b) Investigação 

c) Montagem de instalações 

d) Fabrico de equipamentos 

e) Exploração e manutenção de instalações 
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f) Projecto de instalações 

£) Projecto de equipamentos 
h) Administração de empresas 
|) Comércio de equipamentos. 


Uma sociedade em que existam apenas activi- 
dades de baixa plenitude, por exemplo, comércio, 
trata-se de uma sociedade em que a profissão 
não está enobrecida. Inversamente, se existem 
actividades de elevada plenitude profissional, 
corresponde a um embiente social em que a pro- 
fissão é exercida com grande aceitação. 

Para satisfazer as actividades enunciadas é pos- 
tulado que a sociedade esteja em condições de as 
realizar. Daí a necessidade rígida de evolução 
social: para se permitirem efectuar montagens 
de instalações e respectivos projectos (actividades 
4 e 7) exige-se a laboração de fábricas para o 
fabrico de equipamentos e seus estudos (5 e 8), 
impõe-se a existência de orgãos para exploração 
e manutenção assim como administaação (6 e 9), 
e finalmente é imperativo o intercâmbio de pro- 
dutos pela comercialização dos bens de produ- 
ção (10). 

Uma sociedade em que se verifique este meca- 
nismo é uma sociedade em evolução. Mas o que 
mantém esse mecanismo? Só as actividades de 
grande plenitude profissional poderão impor a 
continuidade da profissão. E surgem assim a inves- 
tigação e a pedagogia na nobre missão de gera- 
dores de plenitude. Para além delas existe o 
transcendente. 
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INSTITUTO NACIONAL DE INVESTIGAÇÃO INDUSTRIAL 


GRUPO DE FORMAÇÃO E APERFEIÇOAMENTO 


Ciclo de Promoção da 


A crescente necessidade de garantir-se um elevado 
nível de rentabilidade num contexto cada vez mais 
competitivo tem obrigado as empresas dos nossos dias 
a dedicarem um cuidado cada vez maior à preparação, 
à coordenação e ao controle de acções que visem o 
desenvolvimento e a adaptação dos seus recursos 
humanos. 

Este ciclo foi preparado para permitir aos promotores 
dessas acções conceber programas eficientes, conduzir 
reuniões de formação segundo técnicas pedagógicas 
adaptadas a adultos e controlar a sua realização. 
Interessa, portanto, aos chefes de pequenas e médias 
empresas, aos directores de pessoal ou serviços de 
formação bem como aos próprios condutores das 
acções, quer sejam chefes operacionais quer monitores 
de formação, 


PFL, 2 
À preparação pedagógica de monitores 
12-16 e 19-23 de Outubro 


das 9.00-1300 h. 
Duração 40 horas 


PFE. 3 


Preparação de monitores de formação 

LISBOA 23-27 de Novembro 
9.30-12.30 € 14.90-17.90 h 
Duração: 30 horas 


MO. 5 
Aplicações estatisticas na actividade da empresa 
Metodologia 


— Cálculo das probabilidades. Teoria dos conjuntos. 
Análise combinatória (revisão). Fenómenos deter- 
minísticos e aleatórios. Teoria das probabilidades, 
Teoria axiomática das probabilidades. 

— Variáveis aleatórias. Definições e propriedades. 
Funções de uma variável aleatória, Caracterização 
de variáveis aleatórias. 

— Distribuições teóricas importantes. Distribuições 
descontínuas. Distribuições contínuas. 

— Estatística descritiva. Apresentação dos dados e sua 
representação gráfica, Determinação dos valores 
característicos mais importantes duma distribuição, 

-— Teoria da amostragem. Definições e conceitos 
básicos, Teorema do limite central. Distribuições 
de funções da amostra: da média, da variância, da 
média de amostras duma população normal, de 
qui-quadrado, de variância empírica, de Student, 
de Snedecor. Controle da qualidade. 

— Ensaios de hipóteses. Introdução, Definições. Ensaios 
de hipóteses simples. Ensaios de hipóteses compos- 
tas. Aplicações do critério de razão de verosi- 
melhanças, 
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Formação na Empresa 


— Estimação. Estimação pontual. Estimação por inter- 
valos. Intervalos de confiança, 

— Correlação linear. 

— Regressão linear. 

— Teoria dos erros de observação. 

Lisboa — 25 — 29 Maio 

10,00 — 13.00 € 15.00 - 18,00 h 

Duração — 30 horas 


MO. 6 
Aplicações estatisticas na actividade da empresa 
Estudos de Casos Reais 


Curso essencialmente prático com a apresentação de 
exemplos diversos dentro dos seguintes temas: 


— Construção de quadros de frequências e correlação. 
Cálculo de medidas descritivas empíricas: médias, 
variância, mediana e moda, 

— Correlação e regressão simples. Erros de previsão. 

— Determinação da dimensão de amostras e estratifi- 
cação de populações. 

— À utilização de conceito de valor esperado na reso- 
lução de problemas de gestão de stoks. 

— Ensaios de hipóteses paramétricas. À importância do 
ensaio do qui-quadrado. 

— Árvores de decisão. 

— Controle estatístico da qualidade. Diagramas de con- 
trole. 

— Sucessões cronológicas. Determinação e eliminação 
da tendência. Identificação e correcção das flutuações 
sazonais. Existência de flutuações cíclicas. 

— Números índices. Sua importância nas aplicações 
correntes e suas limitações. 


, Ciclo de Informática e Gestão 
G. 3 


Introdução à análise de sistemas 
D — Estudos e implantação 


Após um conhecimento considerado suficiente dos 
equipamentos que asseguram a realização do trata- 
mento automático da informação, seguir-se-á em último 
capítulo, a descrição dos princípios que devem infor- 
mar a análise de qualquer sistema no sentido de, 
criando-o em condições de satisfação para as empre- 
sas utilizadoras, permitir tirar o maior partido dos 
equipamentos que o servem. 
Serão apreciadas em pormenor as etapas de um pro- 
jecto de Análise de Sistemas, terminando este ciclo de 
seminários pela descrição de algumas aplicações prá- 
ticas correntes, que servirão para concretizar os conhe- 
cimentos adquiridos. 

D 1 — Estudos de Viabilidade 

Dz — Análise de Detalhe 

D 3 — Lançamento em Execução 
LISBOA — g-11 Novembro 
9 00-12-30 e 14.30-18,00 
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ANTI-CORROSIVA ESPECIAL 


Para a pintura de: 


DEPÓSITOS DE ÁGUA POTÁVEL E OUTROS 
CALDEIRAS 

MÁQUINAS 

TURBINAS 

ESTRUTURAS ABAIXO DOS PAVIMENTOS 
ISOLAMENTOS DAS TUBAGENS 

ESTRUTURAS METÁLICAS 

MATERIAL CIRCULANTE 

GUINDASTES 

POSTES 

PAVIMENTOS DE AÇO, FERRO, MADEIRA, TIJOLO, 
PEDRA E CIMENTO EXPOSTOS AO TEMPO, ETC. 


GRANDE RESISTÊNCIA AO CALOR 


REPRESENTANTE: 


"MENDES DE ALMEIDA 


COMÉRCIO E INDÚSTRIA, S. A. R.1L. 


AV. 24 DE JULHO, 52 LISBOA-PORTUGAL TEL. 667710 


TECNICA XXI 


ak E 


O maior fabricante mundial de rolamentos 
oferece-lhe agora, também outros produtos, 
dos quais destacamos: 


Acoplamentos «KOPPERS» 
Ferramentas de roscar 

Grupos cônicos 

Mangas de esferas «STAR» 

Massa lubrificante 

Material têxtil 

Pontos rotativos 

Redutores velocidade «SINCRON» 
Retentores especiais «V-RING » 
Retentores para indústria «RUBER» 
Rodas e rodizios 

Rolamentos de esfer as; agulhas e roletes 
Rótulas lisas 

Tubos para condução de líquidos corrosivos incluindo água salgada 
Variadores «SINCRON» 


apoiados como sempre numa 
ASSISTÊNCIA TECNICA PERFEITA 


SOCIEDADE E LIMITADA 
Praça da Alegria, 66-A — LISBOA R. Delfim Ferreira, 604 —- PORTO 
Teleirnes 36 2301-2-3-4 Telefones 692054-5-6-7-8 
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